AD-A189  738  CETANE 


uavne  state  univ 

8ESEARCH  N  A  HE; 
AAG29-84-K-8164 


:  SCALE  EVALUATION  AND  POSS  BLE  NODI 
STATE  UNI U  DETROIT  HI  CENT  R  FOR  AU 
:CH  N  A  HENEIN  23  DEC  87  AlO-21397. 


AD- A 189  730 


N.A.  HENEIN 

PROFESSOR  OF  MECHANICAL  ENGINEERING 
DIRECTOR,  CENTER  FOR  AUTOMOTIVE  RESEARCH 
WAYNE  STATE  UNIVERSITY 
DETROIT,  MI 


23  DECEMBER  .1987 
U.S.  ARMY  RESEARCH  OFFICE 
CONTRACT  NO.  DAAG  29-84-K-0164 


WAYNE  STATE  UNIVERSITY 
COLLEGE  OF  ENGINEERING 
CENTER  FOR  AUTOMOTIVE  RESEARCH 
DETROIT,  MI  48202 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


DTIC 


ELECTE 
FEB  021988 


8«  1  *7  O74 


CETANE  SCALE  EVALUATION  AND  POSSIBLE  MODIFICATION 

FINAL  REPORT 


N.A .  HENEIN 

PROFESSOR  OF  MECHANICAL  ENGINEERING 
DIRECTOR,  CENTER  FOR  AUTOMOTIVE  RESEARCH 
WAYNE  STATE  UNIVERSITY 
DETROIT,  MI 

23  DECEMBER  1987 
U.S.  ARMY  RESEARCH  OFFICE 
CONTRACT  NO.  DAAG  29-84-K-0164 


WAYNE  STATE  UNIVERSITY 
COLLEGE  OF  ENGINEERING 
CENTER  FOR  AUTOMOTIVE  RESEARCH 
DETROIT,  MI  48202 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


QU< 

.  iNSP^ca-i) 


i  Accession  For _ . 

~NTIS  GRAM  W 

DTIC  TAB 

Unannouncod 
Juat  1 1'  1  cat  ion - 

By .  . . — 

D1 3tribut Ion/ 
Availability  Codes 


Avail  and/or 


□  □ 


The  view,  opinions,  and/or  findings  contained  in  this  report 
are  those  of  the  author(s)  and  should  not  be  construed  as  an 
official  department  of  the  army  position,  policy,  or  decision, 
unless  so  designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 


2a.  SECURITY  CIASWKA 


2b.  DECULSSIFICATION  /  DOWNGRADING  SCHEDU 


A  PERFORMING  ORGANIZATION  REPORT  NUMBERS) 


3  .  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER($) 


Pinal 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Wayne  State  University 
College  of  Engineering 


6c  ADDRESS  (Oty,  Stott.  and  ZIP  Coda) 

Detroit,  MI  48202 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  appikabla) 

Final  U.  S.  Army  Research  Office 


7b.  ADDRESS  (Oty.  State,  and  HP  Coda) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8a.  NAME  OF  FUNDING  1  SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 


8c  AOORESS  (Oty.  State,  and  HPCoda) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8b.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
Of  applicable) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


H.  TITLE  Onduda  Security  Ctasufkabon) 

CETANE  SCALE  EVALUATION  AND  POSSIBLE  MODIFICATIONS 


12  PERSONAL  AUTHORS) 

Heneln.  Naelm  A‘. 


EFHP 


5.  PAGE  COUNT 

108 


16.  SUPPLEMENTARY  NOTATION 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those 

of  the  Ar®y  position. 


17. _  COSATI  COOES  I  18.  SUBJECT  TERMS  (Continue  on  raaana  if  nataaary  and  idandfy  by  block  numbar ) 

0810  I  G*ouf  I  iu^GROup  I  Cetane  number,  diesel  fulls,  ignition  delay. 


cold  starting. 


(Continue  on  raaam  If  netesury  and  idandfy  by  blodr  numbar) 

An  evaluation  is  made  of  the  validity  and  accuracy  of  the  Cetane  Scale 
in  rating  fuels  nsed  in  diesel  engines.  Six  problems  have  been  identified, 
analyzed  and  possible  solutions  have  been  proposed.  Three  of  the  problems 
are  intrinsic  to  the  scale  and  are  related  to  its  inability  to  rate  low 
ignition  quality  fuels,  instability  of  the  ignition  delay  meter  and  the 
inconformity  of  the  CFR  combustion  system  with  actual  diesel  systems.  The 
other  three  problems  are  related  to  fuels  and  engines.  These  are  the 
invalidity  of  C.N.  in  correlating  the  ignition  delay  data,  its  inaccuracy 
in  ranking  overall  engine  performance  and  emissions  and  its  inability  in 
rating  the  cold  startability  of  engines  on  different  fuels.  The  cold  starta- 
bility  problem  is  investigated  in  two  engines,  under  steady  and  actual 
starting  conditions  in  a  cold  room.  Also,  a  mathematical  model  is  developed 
to  predict  the  startabllitfryof  the  engine.  These  investigations  indicate 


20.  DISTRIBUTION  /  AVAILABILITY  OP  ABSTRACT 
□  UNCLASStfNOAJNUMrrED  □  SAME  AS  RPT. 


22b.  NAME  OP  RESPONSIBLE  INOIVKHJAL 


OTIC  USERS 


22b.  TELEPHONE  Onduda  Ana  Coda)  1 22c  OFFICE  SYMBOL 


DO  FORM  1473, 84  MAR 


83  APR  tdWon  may  b«  uaad  until  Mhbuitbd. 
AN  othbr  MMtioni  bfbfltMOlbtb. 


UNCLASSIFIED 


UNCLASSIFIED 

MUTT  CLAlWCATIOtl  Or  TWO  >55 


ery  clearly  that  the  cetane  ragtag  Is  a  hot  engine,  high  temperature 
test  and  Is  unable  to  rate  the  cold  starfcfcbility  of  diesel  engines. 

Csl  W  JL  n^wV.0  X'l  D  V 

-t 


UNCLASSIFIED 

UCUNITV  CLASSIFICATION  OF  TMII 


TABLE  OF  CONTENTS 


SECTION 

1 

2 

3 

4 

5 

6 


TITLE 

List  of  Figures 
Executive  Summary 

Cetane  Scale:  Function,  Problems  and  Possible 
So1'  yions. 

C.N.  Rating  of  Cold  Startability  of  Diesel 
Engines  on  Different  Fuels. 

Effect  of  Fuel  Volatility  on  Autoignition  and 
Combustion  at  Low  Ambient  Temperatures. 

Mathematical  Model  for  Engine  Startability. 

Conclusions  and  Recommendations. 


PAGE 

i 

1 

3 

59 

82 

92 

106 


I 


LIST  OF  FIGURES 


Figure 

Title 

Page 

Section  2:  Cetane  Scale:  Function,  Problems  and 
Possible  Solutions. 

f 

1 

Relationship  between  CwN.  and  C.R.  in  the  Cetane 
rating . 

7 

2 

Effect  of  increasing  I.D.  to  17°C  C.A.  on  the 
C.N.-C.R.  relationship. 

9 

« 

3 

Truncated  conical  experimental  high  compression 
plug. 

10 

4 

Comparison  between  C.N.-C.R.  relationship  with 

150  F  and  350  F  inlet  air  temperature. 

12 

1 

5 

CJrcle-to-cycle  variation  for  low  ignition  quality 
blends  of  different  cetane  numbers. 

14 

6 

Effect  of  fuel  cetane  number  on  cyule-to-cycle 
variation  in  ignition  delay  at  350°F  and  150°F 
air  temperature. 

15 

i 

7 

General  plan  view  of  CFR  cetane  cylinder  head. 

17 

8 

Installation  of  AVE  12QP  300  CVK  water  cooled 
pressure  transducer. 

18 

1 

9 

Injector  and  pre-chamber  in  the  CFR-Cetane  engine. 

19 

10 

Variation  of  pre-chamber  depth  with  C.R.  in  the 
CFR-cetane  engine. 

21 

l 

11 

Effect  of  needle  opening  pressure  dm  I.D.  in  the 
CFR-cetane  engine. 

22 

12 

Effect  of  I.D.  on  mean  temperature  and  pressure 
during  I .  D„ 

27 

1 

13 

Effect  of  C.N.  on  I.D.  at  different  ih&et  air  temp¬ 
eratures. 

29 

14 

Effect  of  C.N.  on  I.D.  at  different 'mean  temperature 
during  I.D. 

30 

< 

15 

Effect  of  mean  pressure  on  I.D. 

31 

16 

Log  I.D.  plotted  vs  — for  fuels  and  blends  of 
different  C.N.S.  m 

33 

1 

17 

Relationship  between  E  and  C.N.  at  Pm*=400  PSIA. 

34 

i 


18  Global  combustion  parameters  for  diesel  engines.  37 

19  Development  of  combustion  at  low  ambient  temperatures  41 

with  DF  2. 

20  Development  of  combustion  failure  at  low  ambient  tempera-  42 
tures  with  Ref  2  faal  (heavy  fuel). 

21  Development  of  combustion  failure  at  low  ambient  tempera-  43 
tures  with  Ref  1  fuel  (volatile  fuel) . 

Section  3’  C.N.  Rating  of  Cold  Startability  of  Diesel 
Engines  on  Different  Fuels . 

22  Plot  of  angular  velocity  vs.  engine  revolutions.  Fuel:  61 

DF  2,  Oil:  arctic  oil,  ambient  temperature: 20 °C,  cold 

soak  time:  5  hrs . ,  startability  conditions  Started- 
easily. 

23  Plot  of  angular  velocity  vs.  engine  revolutions,  fuel:  62 

DF  2,  oil:  arctic  oil,  ambient  temperature : 15 °C,  Cold 

soak  time:  5  hrs.,  startability  condition:  started 
easily . 

24  Plot  of  angular  velocity  vs.  engine  revolutions,  fuels  63 

DF  2,  oil:  arctic  oil,  ambient  temperature:  10°C, 

cold  soak  time:  12  hrs.,  startability  conditions  started 
easily . 

25  Plot  of  angular  velocity  vs.  engine  revolutions,  fuel:  64 

DF  2,  oil:arctic  oil,  ambient  temperature:  5°C»  cold 

soak  time:  hrs.,  startability  condition:  started 

easily . 

26  Plot  of  angular  velocity  vs.  engine  revolutions,  fuel:  65 

DF  2,  oil:  arctic  oil,  ambient  temperature:  2°C,  cold 

soak  time:  5h  hrs.,  startability  condition:  started 
easily . 

27  Plot  of  angular  velocity  vs.  engine  revolutions.  Fuel:  66 

DF  2,  oil:arctic  oil,  ambient  temperature: -5° cold 

soak  time:  13  hrs.,  startability  condition:  started 
easily . 

28  Plot  of  angular  velocity  vs.  engine  revolutions.  Fuel:  67 

DF  2,  oil:arctic  oil,  ambient  temperature:  -8°C,  cold 

soak  time:  12  hrs.,  startability  condition:  with 
difficulty . 

29  Plot  of  angular  velocity  vs.  engine  revolutions,  fuel:  68 

DF  2,  oil:arctic  oil,  ambient  temperature:-10°C,  cold 

soak  time: 5  hrs.,  startability  condition:  did  not  start. 


il 


69 


30  Plot  of  angular  velocity  vs.  crank  angle  ,  fuel:  JP  4 
Bil: arctic  oil,  cold  soak  time:  12  hr»s,  ambient  tempe¬ 
rature:  25°C,  startabillty  condition:  started  easily. 

31  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4,  69 

oil:  arctic  oil,  cold  soak  time:  6  hrs.,  ambient 
temperature: 10 °C,  startabillty  condition:  started  easily. 

32  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4  70 

oil:  arctic  oil,  cold  soak  time:5hrs.,  ambient  tempera¬ 
ture:  5#C,  startabillty  condition:  started! easily  . 

33  Plot  of  angular  velocity  vs.  crank  angle  fuel,  JP4  71 

oil: arctic  oil:  cold  soak  time:  12  hrs . ,  ambient  tempe¬ 
rature:  0°C,  startabillty  condition:  started  easily. 

34  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4,  72 

oil: arctic  oil,  cold  soak  time:  5h  hrs.,  ambient  tempe¬ 
rature  :-5*C,  startabillty  condition:  started  easily. 

35  Plot  of  angular  velocity  vs.  Crank  angle,  fuel:  JP4,  73 

oil:  arctic  oil,  cold  soak  time: 6  hrs.,  ambient  tempe¬ 
rature:^*^,  startabillty  condition:  started  easily. 

36  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4,  74 

oil: arctic  oil,  cold  soak  time:  6  hrs.,  ambient  tempe¬ 
rature  :-8°C,  startabillty  condition:  started  easily. 

37  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4  75 

oil: arctic  oil,  cold  soak  time : S^hrs . ,  ambient  tempe- 
rature:-9°C,  startabillty  condition:  started  easily. 

38  Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4,  76 

oil:  arctic  oil,  cold  soak  time:  73shrs . ,  ambient  tempe¬ 
rature  :-10°C ,  startabillty  conditions:  started  with 
difficulty . 

39  Plot  of  angular  velocity  vs.  crank  angle.  Fuel:  JP4,  77 

oil:  arctic  oil,  cold  soak  time:  6  hrs.,  ambient  tempe¬ 
rature:  -11*C,  startabillty  condition:  started  with 
difficulty . 

40  Flot  of  angular  velocity  vs.  crank  angle.  Fuel:  JP4,  78 

oil:  arctic  oil,  cold  soak  time:8hrs.,  ambient  tempe¬ 
rature:  -13°C,  startabillty  condition:,  failure  to  start. 

41  Comparison  between  the  engine  revolutions  to  maaimitm  80 

starting  speed  on  DF  2  and  JP  4  fuels . 

Section  4 :  Effect  of  Fuel  Volatility  on  Autoignition 
and  Combustion  at  Low  Ambient  Temperatures. 

42  Layout  of  experimental  setup.  84 

43  Schematic  of  electronic  instrumentation.  86 


iii 


44  Variation  of  misfiring  temperature  with  cetane  number  88 

for  distillated  and  blend  fuels. 

45  Variation  of  ignition  delay  period  with  inlet  air  tem-  89 

perature  at  Pm=400  psia  (2.76  MPa),  for  distillated 

and  blend  fuels. 

46  Variation  of  ignition  delay  period  with  cetane  number  90 

at  different  ambient  temperature  (distillated  and  blend 
fuels) . 


iv 


Jk 


SECTION  1 


EXECUTIVE  SUMMARY 


Many  problems  are  experienced  by  the  military  and  industry 
in  using  the  Cetane  Scale  to  rate  the  fuels  for  diesel  engines . 

The  main  objectives  of  this  program  are: 

1.  To  evaluate  the  validity  and  accuracy  of  the  Cetane 
Scale  in  rating  fuels  used  in  diesel  engines. 

2.  Identify  the  problems. 

3.  Propose  solutions. 

4.  Investigate  the  C.N.  as  a  cold  startability  in¬ 
dicator. 

To  meet  the  first  three  objectives  the  function  of  the 
Cetane  Scale  is  identified  and  six  problems  are  examined  and 
solutions  proposed.  These  problems  are: 

1.  Inability  of  the  ASTM  D-613  Method  to  Rate  Low  Ignition 
Quality  Fuels. 

2.  Instability  of  the  Ignition  Delay  Meter  in  Rating  Low 
Ignition  Quality  Fuels. 

3.  Inconformity  of  the  CFR  Combustion  System  with  actual 
Diesel  Systems.  . 

4.  Invalidity  of  C.N".  in  Correlating  the  I.D.  for  Different 
Fuels,  Engines  and  Operating  Conditions. 

5.  Inaccuracy  of  C.N.  in  Ranking  the  Overall  Performance  and 
Emission  Behaviour  of  Fuels  in  Diesel  Engines. 

6.  Inability  of  C.N.  in  Rating  the  Cold  Startability  of 
Diesel  Engines  on  Different  fuels. 

A  possible  solutions  for  these  problems  are  proposed  in  this  report. 

To  meet  the  fourth  objective  experimental  investigations  sup¬ 
ported  by  a  mathematical  model  have  been  conducted.  Experiments 
were  conducted  at  low  ambient  temperatures  in  a  cold  room  at 
steady  operating  conditions  and  under  actual  starting  conditions. 
The  tests  covered  ambient  temperatures  ranging  from  25°C  to  -  45°C 
or  till  the  engine  failed  to  fire  or  start. 
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The  steady  state  experiements  were  conducted  on  a  single 
cylinder  TACOM  engine  on  a  variety  of  fuels.  The  ignition 
delay  and  the  misfiring  temperature  were  measured  for  each 
fuel.  The  results  showed  that  the  ignition  delay  of  the  high 
volatility  fuels  is  shorter  than  the  ignition  delay  for  the 
low  volatility  fuels  having  the  same  Cetane  number.  Also, 
the  misfiring  temperature  was  found  to  decrease  as  the 
volatility  of  the  fuel  increases.  These  results  indicate 
that  the  fuel  volatility  plays  a  major  role  in  the  combustion 
process  of  the  diesel  engine  at  low  ambient  temperatures. 

The  starting  tests  were  conducted  on  air  cooled,  direct 
injection  diesel  engine  in  the  cold  room.  The  results  in¬ 
dicated  that  the  Cetane  number  is  not  capable  of  rating  the 
startability  of  the  engine  on  different  fuels.  A  low  Cetane 
fuel  which  is  volatile  was  able  to  start  the  engine  while 
a  high  Cetane  less  volatile  fuel  failed  to  start  the  engine. 

The  above  investigations  indicate  very  clearly  that  the 
C.N.  is  a  hot  engine,  high  temperature  test  and  is  unable  to 
rate  the  cold  startability  of  diesel  engines. 

It  is  recommended  that  an  investigation  be  made  to  find 
the  necessary  modifications  in  the  scale  to  make  it  of  value 
in  the  rating  of  the  cold  starting  of  the  diesel  engines  on 
different  fuels. 
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SECTION  2 


CETANE  SCALE:  FUNCTION*  BROBLEMS  AND  POSSIBLE  SOLUTIONS 
2.1  BACKGROUND 

Many  efforts  have  been  made  during  the  last  60  years  to  rate  the 
behaviour  of  different  fuels  in  diesel  engines.  Combustion  noise  was 
one  of  the  performance  parameters,  of  concern  to  researchers.  As  early 
as  1931  LeMesurier  and  Stansfield  [1]  proved,  from  experiments  on 
several  compression-ignition  engines,  that  the  combustion  knock  was 
broadly  related  to  the  delay  which  occurred  between  the  moment  of  injec¬ 
tion  and  the  beginning  of  rapid  pressure  rise.  They  reported  that 
different  engines  seemed  to  have  only  slightly  different  relative 
effects  on  the  fuels  tested.  Subsequent  investigators  such  as  Boerlage 
and  Broeze  1931  (2] ,  Pope  and  Murdock  1932  [3] ,  Schweitzer  et  al  1933- 
1937  [4,5,6]  used  the  ignition  delay  as  the  criteria  to  rate  the  be¬ 
haviour  of  the  different  fuels  in  diesel  engines. 

The  reference  fuel  of  the  highest  rating  of  100,  on  the  scale  has 
been  n-hexadecane ,  Rbowb  as  n-Cetane.  The  reference  fuel  of  the 
lowest  rating  of  zero  on  the  scale  had  been  alpha-methyl-naphthalene, 
C11H10’  from  1932  t0  1962 '  The  C'N'  accordill8  to  this  scale  is 

C.N.  =  %  n-Cetane  (1) 

In  1962  ASTM  [7]  replaced  alpha-methyl-naphthalene  (<*MN)  with  a  new 
low-cetane  number  reference  fuel,  hepta  methylnonane  (HMN).  The  cetane 
rating  for  heptamethylnonane  is  15 . 

The  definicion  of  the  ASTM  C.N.  [7]  is  the  whole  nomber  nearest  to 
the  value  determined  by  calculation  from  the  percentage  by  volume  of 
normal  cetane  in  a  blend  with  HMN  which  matches  the  ignition  quality  of 
the  test  fuel  when  compared  by  the  D-613  ASTM  Cetane  Method.  The  match- 
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ing  blend  percentages  to  the  first  decimal  are  inserted  in  the  following 
equation  to  obtain  the  C.N. 


C.N.  =  X  n-cetane  +  0.15  (%  HMN)  (2) 

The  single  cylinder  engine  used  in  the  fuel  rating  is  known  as  CFR 
(Cooperative  Fuels  Research)  engine,  manufactured  by  Waukesha  Company. 
The  engine  specifications  are  given  in  Appendix  A,  and  the  operating 
conditions  for  the  ASTM  D-613  test  procedure  are  given  in  Appendix  B. 

There  has  been  growing  concern  during  the  last  two  decade  over  many 
problems  associated  with  the  Cetane  number  determined  by  the  ASTM  D-613 
method,  particularly  for  alternate  fuels  which  have  been  considered  as 
replacements  of  the  regular  diesel  fuel  or  as  fuel  extenders  [8]. 
Additional  concern  is  over  the  use  of  C.N.  to  rate  the  performance, 
gaseous  emmissions,  exhaust  solid  particulates,  white  smoke  and  the  cold 
startability  of  diesel  engines. 

The  purpose  of  this  section  is  to  identify  the  ortRin.il  function  of 
the  Cetane  rating  of  the  fuel  by  using  the  ASTM  D-613  method,  and  to 
examine  the  different  problems  of  concern  to  the  technical  community. 

This  section  examines  six  problems  and  introduces  nossible  solutions 

or  recommends  actions  to  reach  solutions.  These  problems  are: 

1.  Inability  of  the  ASTM  D-613  Method  to  Rate  Low  Ignition  Qua¬ 
lity  Fuels. 

2.  Instability  of  the  Ignition  Delay  Meter  in  Rating  Low  Ignition 
Quality  Fuels. 

Inconformity  of  the  CFR  Combution  System  with  Actual  Diesel 
Systems . 
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4.  Invalidity  of  C.N.  in  Correlations  the  I.D.  for  Different 


Fuels,  Engines  and  Operating  Conditions. 

5.  Inaccuracy  of  C.N.  in  Ranking  the  Overall  Performance  and 
Emission  Behaviour  of  Fuels  in  Diesel  Engines. 

6.  Inability  of  C.N.  in  Rationg  the  Cold  Startability  of  Diesel 
Engine  on  Different  Fuels. 

2.2  FUNCTION  OF  C.N. 

The  basic  function  of  the  C.N.  is  to  rate  the  autoignition  proper¬ 
ties  of  fuels  in  diesel  engines.  The  details  of  the  autoignition 
chemical  processes  are  not  known  because  they  are  very  complex,  fast  and 
involve  many  unstable  combustion  radicals  which  are  not  easy  to  trace. 
The  autoignition  chemical  processes  are  preceeded  by  physical  processes 
which  are  dependent  on  the  fuel  physical  properties  and  the  design  of 
the  injection  system  and  the  combustion  chamber  of  the  engine.  In  the 
ASTM  D-613  the  measurement  of  the  autoignition  delay  of  the  fuel  is  done 
by  measuring  the  ignition  delay  in  the  standard  engine  under  standard 
operating  conditions,  and  adjusting  the  compression  ratio  C.R.  such  that 
the  I.D.  is  2.41  m  sec.,  or  13°  CA  at  900  RPM.  The  start  of  injection 
is  at  13°  CA  before  TDC  and  the  pressure  rise  due  to  combustion  (or  the 
end  of  I.D.)  should  take  place  at  TDC. 

No  other  combustion  or  engine  performance  parameter  is  measured 
during  the  rating  test. 

2.3  PROBLEMS  AND  POSSIBLE  SOLUTIONS 

In  this  section  six  problems  are  examined 

The  first  three  are  intrinsic  problems  associated  with  ASTM  D-613 
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procedure  or  the  CFR  engine  design.  The  last  three  are  fuel/engine 
probleas  related  to  the  behaviour  of  fuels  in  actual  diesel  engine. 


\ 

2.3.1  Problem  1:  Inability  of  ASTM  D-613  Method  to  Rate  Low  Ignition 
Quality  Fuels: 

Very  high  compression  ratios  are  needed  to  rate  alternate  fuels, 
haveing  low  cetane  numbers.  Figure  1  shows  the  relationship  between  the 
compression  ratio  and  C.N.  according  the  ASTM  D-613  and  equation  (2)  in 
which  HMN  is  used  as  the  bottom  of  the  scale  at  15  C.N.  The  same  figure 
shows  the  same  relationship  with  aMN  as  the  zero  C.N.  fuel  according  to 
equation  (1).  It  is  clear  that  the  relationship  is  close  to  linear  for 
fuels  haveing  C.N.  between  100  and  50.  At  C.N.  less  than  50  the  slope 
of  the  curve  increases  at  a  high  rate  as  the  C.N.  decreases.  Figure  1 
shows  that  a  C.R.  of  about  29.5  is  needed  for  rationg  a  fuel  having  C.N. 
=  10.  Higher  compression  ratios  are  needed  for  rating  lower  C.N.  fuels. 

The  very  high  compression  ratios  needed  for  rating  low  C.N.  fuels 
has  been  the  concern  of  many  investigators  [9].  Urban  and  Gray  1968 

[10]  of  the  U.S.  Army  Fuels  and  Lubricants  Research  Laboratory  reported 
the  results  of  a  comprehensive  study  on  the  feasiblity  of  three  alter¬ 
nate  methods  to  solve  this  problem  and  rate  low  ignition  quality  fuels. 
The  first  method  was  to  blend  n-cetane  fuel  with  the  sample  and  then 
calculate  the  sample  C.N.  from  the  C.N.  of  the  blend.  This  method  was 
found  to  be  unsatisfactory.  More  recently,  Hardenberg  and  Ehnert  1981 

[11]  reported  dissatisfaction  with  this  method  in  determining  the  C.N. 
of  alcohols.  The  second  method  reported  by  Urban  and  Gray  [10]  is  to 
correlate  the  C.N.  and  Octane  number  of  the  sample.  They  reported  that 
the  C.N.  estimated  by  this  method  was  less  raliable  than  C.N.  directly 


determined.  The  third  method  is  to  modify  the  ASTM-D613  procedure  by 
using  longer  I.D.  than  the  standard  delay  of  13°  C.A.  at  900  RPM. 

Figure  2  shows  the  results  obtained  for  the  present  program  and  indicate 

\ 

a  small  reduction  in  the  C.R.  as  a  result  of  using  an  I.D.  of  17°  C.A. 
instead  of  the  standard  13°  C.A.  Similar  results  were  reproted  by 
Bowder  [12]. 

Another  method  to  solve  this  problem,  proposed  by  the  French  Army 
[13]  to  Waukesha  Motor  Company,  is  to  use  a  compression  plug  which  has  a 
notch  cut  deeply  into  one  side  to  clear  the  transfer  port  between  the 
two  chambers  in  order  to  increase  the  highest  C.R.  of  the  standard  CFR 
engine.  This  proposal  was  not  favored  by  Waukesha  Motor  Company  because 
of  the  tendency  of  this  plug  to  create  a  different  flow  pattern  and  tur¬ 
bulence  in  the  prechamber  of  the  standard  CFR  engine  and  the  different 
machining  operations  required  for  the  tapered  shape  of  the  slot. 
Waukesha  proposed  back  to  the  French  Any  the  use  of  a  modified  high 
compression  ratio  plug  which  has  a  symmetrical  truncated  conical  end,  to 
reach  C.R.  higher  than  31:1.  Figure  3  shows  such  a  plug.  No  tests  were 
reported  on  the  effectiveness  of  this  plug  in  rating  the  low  ignition 
quality  fuels. 

An  analysis  of  the  parameters  which  affect  the  autoignition  process 
of  the  fuels  in  any  diesel  engine  [14]  indicated  that  the  I.D.  is  most 
sensitive  to  the  air  teaiperature  and  molecular  structure  of  the  fuel,  in 
addition  to  the  air  pressure  according  to 
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(3):  Truncated  conical  experimental  high 
compression  plug. 


T 


where  E  is  the  Global  activation  energy  for  the  preignition  reactions 
and  depends  mainly  on  the  molecular  structure  of  the  fuel. 

Both  the  air  temperature  and  pressure  increase  with  C.R.  in  the  CFR 

\ 

engine.  The  relative  effects  of  temperature  and  pressure  on  the  igni¬ 
tion  delay  are  discussed  in  details  in  a  previous  publication  [14].  The 
effect  of  the  increase  in  temperature  with  C.R.  is  much  more  than  the 
effect  of  the  increase  in  pressure  on  the  l.D.  This  has  been  found  to 
be  particularly  true  for  low  C.N.  fuels  which  have  relatively  higher 
activation  energies.  Based  on  this  analysis,  it  becomes  clear  that 
raising  the  inlet  air  temperature  above  the  150°F,  used  in  the  0-613 
method,  would  reduce  the  C.R.  needed  to  produce  an  l.D.  of  2.41  m  sec. 
for  the  low  ignition  quality  fuels.  This  has  been  verified  experi¬ 
mentally  by  Elias  [15].  Figure  4  shows  2  cetane  scales:  the  original 
cetane  scale  at  determined  by  using  ASTM  D-613  method,  with  an  air  inlet 
temperature  of  150°F,  and  the  scale  developed  after  increasing  the  inlet 
air  temperature  from  150°F  to  350°F.  With  the  new  scale  it  is  possible 
to  rate  low  ignition  quality  fuels  without  reaching  the  very  high  C.R.s 
needed  in  the  original  scale.  A  zero  C.N.  fuel  can  be  autoignited  with 
an  l.D.  of  2.41  m  sec.,  at  900  RPM  at  a  C.R.  =  19.4:1. 

The  solution  to  this  problem  is  to  raise  the  inlet  temperature 
above  the  150°F  specified  by  the  ASTM  0-613  procedure.  This  requires 
the  use  of  an  electric  air  heater  of  a  higher  capacity  than  the  heater 
currently  used  on  the  CFR  engine. 

Ratings  with  inlet  air  temperatures  between  150°F  and  350°F  need  to 
done  to  choose  a  temperature  which  can  be  reached  by  using  an  electric 
heater  of  a  reasonable  capacity. 
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Pig.  (4):  Comparison  between  C.N.-C.R.  rela¬ 
tionship  with  150®F  and  350°F  inlet 
air  temperatures. 
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2.5.2 


Problem  2:  Instability  of  the  Ignition  Delay  Meter  in  Rating 


Low  Ignition  Quality  Fuels: 

Many  investigators  observed  that  the  ignition  delay  meter  becomes 

\ 

unstable  and  performs  erratically  as  the  C.N.  of  the  sample  drops  below 
twenty.  Urban  and  Gray  (10]  proposed  the  use  of  an  oscilloscope  and 
associated  pressure  transducers  instead  of  the  standard  ASTM  ignition 
delay  meter,  to  measure  the  I.D.  of  the  low  C.N.  fuels.  They  reported 
that  the  accuracy  of  the  rating  depended  on  the  operator's  ability  to 
average  continuous  cycles  on  the  oscilloscope  due  to  cycle-to-cycle 
variations  in  the  pressure  trace. 

The  reason  for  such  instability  has  not  been  known  till  Sczomak 
(16]  proved  experimentally  that  it  is  caused  by  the  cycle-to-cycle 
variation  in  the  autoignition  of  low  C.N.  fuels.  Figure  5  shows  traces 
for  the  gas  pressure  in  the  prechamber,  the  rate  of  change  of  the  gas 
pressure,  the  needle  displacement,  the  needle  velocity  and  the  fuel 
injection  pressure  before  the  injector.  Each  of  the  figures  (Sa)  to 
(5f)  is  for  twenty  superimposed  consecutive  cycles.  It  is  noticed  that 
the  variations  in  the  cylinder  pressure  are  caused  by  variations  in  the 
I.D.,  and  these  variations  increase  as  the  C.N.  decreases.  It  should  be 
noted  that  the  other  traces  for  the  needle  displacement,  needle  velocity 
and  fuel  injection  pressure  showed  no  cycle-to-cycle  variations  for  all 
the  fuels. 

Many  attempts  have  been  made  to  reduce  the  cycle-to-cycle  varia¬ 
tions  in  I.D.,  which  result  in  the  instability  of  the  I.D.  meter. 
Figure  6  from  reference  (16]  shows  the  effect  of  increasing  the  inlet 
air  temperature  from  150°F  to  350AF  on  the  coefficient  of  variation  in 
the  I.D.  It  is  clear  that  the  cycle-to-cycle  variations  are  reduced  for 
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C.N.s  less  than  18.  However,  these  variations  still  exist  even  at  this 
very  high  inlet  air  temperature. 

At  the  present  time  the  processes  which  result  in  the  cycle-to- 
cycle  variations  in  the  autoignition  process  are  not  known.  Further 
investigations  in  this  area  are  needed. 

2.3.3  Problem  3:  Inconformity  of  the  CFR  Combustion  System  with 

Actual  Diesel  Combustion  Systems: 

The  autoignition  of  fuel  sprays  injected  in  an  engine  is  the  result 
of  two  main  sets  of  processes:  physical  and  chemical.  The  physical 
processes  lead  to  the  formation  of  a  combustible  mixture,  with  an  equi¬ 
valence  ratio  close  to  the  stoichiometric  ratio.  The  chemical  processes 
are  mainly  decomposition  processes  leading  to  the  formation  of  the 
combustion  redicals  and  the  oxidation  processes  which  release  the  energy 
of  the  combustion  reactions:  In  actual  diesel  engines  the  combustion 
systems  are  designed  to  allow  the  above  processes  to  take  place  and  keep 
the  ignition  delay  within  certain  limits  such  that  the  majority  of  the 
combustion  energy  is  released  close  to  TDC  of  the  engine.  This  requires 
the  optimization  of  the  combustion  system,  which  consists  of  the  injec¬ 
tion  system  and  combustion  chambers,  to  suit  the  properties  of  the  fuel. 
The  injection  parameters  include  the  nozzle  type  and  size,  injection 
pressure,  and  injection  timing. 

In  the  CFR-Cetane  engine  the  shape  of  the  combustion  chamber  varies 
according  to  the  C.R.  which  in  turn  varies  according  to  the  ignition 
quality  of  the  fuel  being  rated.  Figures  7,8,  and  9  give  different 
views  of  the  CFR  prechamber,  main  chamber  and  injector.  The  C.R.  of  the 
engine  is  controlled  by  varying  the  volume  of  the  prechamber  by  moving 
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FIG.  8 

CROSS- SECTION  B-B  OF 
FIGURE  (  6  ).  SHOWING 
INSTALLATION  OFAVLI2QP 
300CVK  WATEP  COOLI'D 
PRESSURE  TRANSDUCERS 
IN  BOTH  PRE  AND  MAIN 
CHAMBERS-  (25) 
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(9):  Cross-section  A-A  of  figure  6  showing  the 
injector  and  pre-chamber  in  the  CFR-Cetane 
engine . 
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the  compression  plug  in  the  prechamber.  The  distance  between  the  flat 
surface  of  the  compression  plug  and  injector  nozzle  tip  decreases  from 

2.435"  at  C.R.  =  7.92  to  0.435"  at  C.R.  =  31.  Figure  10  from  reference 

\ 

[15]  shows  the  shape  of  the  prechamber  at  different  C.R.s.  An  estima¬ 
tion  of  the  spray  penetration  indicated  that  the  sprays  impinges  on  the 
surface  of  the  compression  plug,  at  all  compression  ratios  during  the 
ignition  delay  period  of  2.41  m  sec.  [15].  It  is  not  clear  at  this  time 
whether  ignition  starts  in  the  fuel  spray  before  or  after  it  impinges  on 
the  compression  plug,  particularly  at  the  high  compression  ratios. 

In  actual  diesel  engines  the  ignition  delay  is  mainly  controlled  by 
the  chemical  preignition  processes  rather  by  the  physical  processes  even 
at  the  low  ambient  temperatures  [17].  The  physical  processes  are  mainly 
dependent  on  the  injection  system  and  the  extent  of  spray  atomization. 
For  a  fixed  nozzle  configuration  the  atomization  process  depends  to  a 
large  extent  on  the  injection  pressure,  which  is  controlled  by  the 
needle  opening  pressure.  If  an  increase  in  the  injection  pressure 
produces  a  decrease  in  I.D.,  this  is  an  indication  that  the  atomization 
process  is  not  appropriate. 

Figure  11  shows  the  effect  of  increasing  the  needle  opening  pres¬ 
sure  on  the  I.D.  of  different  fuels  in  the  CFR  engine  at  C.R.  =  13.6 
[20].  The  fuels  used  range  from  light  fuels  n-C^,  typical  Japanese  fuel 
JPN  with  50%  boiling  point  of  263. 5°C,  heavy  fuel  1  with  50%  boiling 
point  of  398°C  and  a  heavier  fuel  2  with  5U%  boiling  point  of  459°C. 
Figure  11  shows  that  the  I.D.  of  the  two  light  fuels  decrease  as  the 
injector  opening  pressure  increases  from  1000  psi  to  1500  psi,  and 
remains  constant  at  higher  needle  opening  pressure.  The  I.D.  for  the 
two  heavier  fuels  decreased  continuously  with  the  increase  in  injection 
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Pig.  (11):  Effect  of  needle  opening  pressure  on 
I.D.  In  the  CPR-Cetane  engine. 
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pressure.  This  indicates  that  part  of  the  increase  in  I.D.  is  due  to 
the  improper  atoaization  process.  The  saae  effect  has  been  reported  by 
Lyn  et  al.,  (18]  and  Pederson  et  al.,  [19]. 

A  possible  solution  of  this  problem  is  to  increase  the  needle 
opening  pressure  for  the  heavier  fuels.  This  would  require  rating  the 
fuels  in  the  CFR  engine  at  different  needle  opening  pressures.  The 
rating  would  be  done  with  a  needle  opening  pressure  beyond  which  the 
rating  is  independent  of  the  fuel  pressure. 

2.3.4  Problea  4:  Invalidity  of  C.N.  in  Correlating  the  I.D.  for 

Different  Fuels,  Engines  and  Operating  Conditions 

There  has  been  concern  over  the  validity  of  the  C.N.  in  rating  the 
ignition  delay  of  fuels  in  actual  diesel  engines.  Many  experiments  have 
been  conducted  recently  in  engines  [21,22,23]  and  constant  volume  bombs 
[24,25]  to  determine  the  I.D.  of  a  wide  range  of  fuels  and  correlate 
I.D.  with  the  C.N  as  determined  by  the  D-613  procedure. 

Doyle  and  Needham  [21,22]  reported  the  results  of  extensive  tests 
on  a  single  cylinder  engine  using  ten  fuels  in  both  IDI  and  DI  combus¬ 
tion  ays  teas .  The  main  fuel  properties  are  listed  in  Appendix  C.  The 
fuels  included  DF-2,  2  aiddle  distillates  of  different  viscosities  and 
aromatic  contents,  a  high  volatility  broad  cut  fuel  composed  of  66% 
Naphtha  and  34%  DF-2,  degummed  sunflower  oil,  a  more  volatile  less 
viscous  sunflower  aonoester,  shale  derived  f-*el,  a  blend  of  58%  DF-2  and 
42%  coal  derived  SRC  II,  a  blend  of  78%  DF-2  and  10%  Methanol  and  12% 
Hexan-1-01,  and  1990  Canadian  Diesel  Fuel  from  tar  sands.  The  fuels 
covered  the  following  ranges:  C.N.:  29  to  50.6,  molecular  weight:  172 
to  635,  Aromatic  contents  of  hydrocarbon  fuels:  19%  to  57%  and  kine- 
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matic  viscosity  at  40°C:  0.91  to  33.93.  The  tests  covered  a  wide  range 
of  engine  operation  conditions:  low  speed,  peak  torque  and  rated 
speeds,  full  load,  light  load  and  idle.  Doyle  and  Needham  reported  that 
there  is  overall  reasonable  agreement  between  I.D.  and  C.N.  in  both  D.I. 
and  I.D. I.  engines  for  petrolium  and  synthetic  derived  mid-distillates 
and  the  sunflower  ester.  However,  I.D.  of  sunflower  oil  and  volatile 
blends  of  Naphtha  and  Methanol  with  DF-2  showed  varying  delay  response 
to  load  and  speed. 

The  sunflower  oil  produced  shorter  I.D.  than  the  other  fuels  having 
the  same  C.N.  at  high  speeds  and  loads.  Two  factors  might  have  contri¬ 
buted  to  this  descrepency.  The  first  in  the  Cetane  rating  of  29  for  the 
sunflower  oil  is  done  in  the  CFR  at  the  standard  D-613  procedure.  As 
indicated  earlier  this  procedure  calls  for  a  constant  injector  opening 
pressure  of  1500  psi  which  is  enough  to  produce  a  proper  spray  for  light 
or  medium  petroleum  derived  fuels.  The  author  feels  that  the  cetane 
rating  of  the  Sunflower  oil  does  not  represent  its  autoignition  quality, 
because  of  the  CFR  injection  system  being  unsuitable  for  sunflower  oil. 
The  sunflower  oil  has  a  molecular  weight  of  635  as  compared  to  187  for 
DF-2.  The  molecular  weight  for  the  petroleum,  shale  or  coal  derived 
fuel  tested  ranges  between  172  to  208.  The  viscosity  of  the  sunflower 
oil  is  33.93  CST,  at  40°C,  while  the  viscosity  of  all  other  fuels  ranges 
between  0.91  and  4.37.  It  is  clear  that  both  the  molecular  weight  and 
viscosity  of  sunflower  oil  are  quite  different  from  the  rest  of  the 
fuels.  Accordingly  it  appears  that  the  C.N.  rating  of  sunflower  oil 
should  take  this  into  consideration  and  adjust  the  injection  parameters 
in  such  a  way  as  to  develop  a  proper  fuel  suray  in  the  CFR  engine.  The 
second  factor  is  related  to  the  optimization  of  the  injection  systems  on 
the  engines  for  the  use  of  sunflower  oil. 
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Siebers  [24]  used  s  constant  voluae  boab  to  measure  the  I.D.  of 
several  fuels,  some  of  which  were  tested  by  Needham  and  Doyle  [22],  as 
explained  earlier.  Other  fuels  were  blends  of  n-Cetane  and  HKN  with 
C.N.  ranging  froa  15  to  100.  A  preaixed  charge  of  02^/112/02/^  mixture 
was  ignited  and  burned  to  fora  the  high  pressure,  high- temperature 
environaent  and  sianilate  the  diesel  theraodynaaic  conditions.  The 
pressure  ranged  froa  10  ata.  to  50  ata. ,  and  the  temperature  ranged  froa 
600°K  to  1400°K.  Siebers  [24]  concluded  that  the  C.N.  cannot  provide  a 
consistent  and  accurate  aeasure  of  I.D.  of  fuels  whose  I.D.  dependence 
on  teaqperature  and  type  of  ignition  process  (single-stage  or  two-stage) 
different  from  those  of  the  reference  fuel  blends  of  n-cetane  and  HMN. 
Sieber  observed  a  pronounced  two  stage  ignition  for  HMN  at  temperatures 
less  than  900°K  at  a  pressure  of  30  atm. 

The  two-stage  ignition  has  been  reported  in  a  recent  study  [20]  to 
take  place  for  fuels  having  C.N.  between  19.4  and  73.  The  low  C.N. 
fuels  produced  the  two  stage  ignition  at  higher  compression  ratios  than 
the  high  C.N.  fuels. 

In  the  author's  opinion,  C.N.  cannot  correlate  the  I.D.  in  all 
types  of  engines,  for  all  fuels  derived  from  different  sources,  under 
all  engine  operating  speed  and  load  conditions.  This  may  be  due  to  the 
fact  that  the  Cetane  rating  is  conducted  only  under  a  set  of  specified 
conditions  in  which  the  I.D.  is  kept  constant  and  the  air  pressure  and 
temperature  are  varied  simultaneously  by  adjusting  the  C.R.  at  a  con¬ 
stant  speed  and  constant  inlet  air  and  coolant  temperatures.  In  actual 
diesel  engines  the  C.R.  is  kept  constant  while  the  ambient  temperature 
varies  for  cold  start  to  hot  start  conditions,  the  intake  air  pressure 
and  temperature  vary  according  to  whether  the  engine  in  naturally 
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aspirated  or  turbocharged,  the  engine  speed  varies  fora  idling  to  rated 
speed  and  the  load  or  F/A  ratio  vary  froa  no  load  to  full  load.  The 
autoignition  process  as  reported  by  Siebers  [24],  even  for  the  same 
fuel,  Bay  vary  froa  a  single  stage  under  some  of  the  above  conditions  to 
a  double  stage  under  other  conditions.  Accordingly  it  seems  that  using 
a  C.N.  value  to  rate  the  I.D.  under  all  the  above  operating  conditions 
is  very  optiaistic. 

The  author  feels  than  the  I.D.  can  be  correlated  according  to  equation 
(3),  in  terms  of  three  basic  parameters: 

i.  Air  temperature, 

ii.  Air  pressure. 

iii.  A  global  activation  energy  for  the  preignition 
reactions. 

The  tesiperatures  and  pressures  in  equation  (3)  should  represent  the 
effective  values  for  the  preignition  cheaical  reactions  which  control 
the  preignition  processes.  The  temperature  at  the  start  of  fuel 
injection  is  often  used  for  such  correlations,  while  it  does  not 
represent  the  teaq>erature  which  affect  the  rate  of  the  chemical  reac¬ 
tions.  Figure  12  is  an  illustration  on  the  variation  in  the  temperature 
during  the  I.D.  of  three  fuels. 


P"(C)  <  *(»)  <  P"(B) 
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and  Tm  and  Pm  depend  on  both  the  injection  timing  and  the  length  of  I.D. 
period  in  crank  angle  degrees. 

The  results  of  tests  conducted  on  a  single  cylinder  research  diesel 
engine  under  controlled  temperatures  and  pressures  during  the  I.D.  are 
shown  in  figure  13.  This  figure  shows  the  effect  of  C.N.  on  the  I.D.  of 
several  distillates  and  blends  of  primary  reference  fuels  in  a  D.I. 
single  cylinder  diesel  engine.  The  data  show  the  results  obtained  at 
inlet  air  temperatures  varying  from  -4°F  (-20°C)  to  350°F  (176. 7°C). 
The  inlet  air  pressure  was  adjusted  to  keep  the  mean  pressure  during  the 
ignition  delay  at  a  constant  value  of  400  psia  (  KPa).  For  this, 
the  inlet  pressure  had  to  be  increased  above  atmospheric  for  tempera¬ 
tures  above  the  normal  room  temperature  and  had  to  be  decreased  below 
atmospheric  for  tests  at  low  ambient  temperatures.  The  fuels  used 
consisted  of  blends  of  primary  reference  fuels  n-cetane  (C.N.  =  100)  and 
aMN  (C.N.  =  0)  and  four  distillates:  DF-2,  J.P4,  Ref.  1  and  Ref.  2 
fuels.  The  properties  of  these  fuels  are  given  in  Appendix  D.  Figure/4 
shows  the  effect  of  changing  the  mean  temperature  during  the  I.D.  for 
fuels  having  different  C.N.s. 

In  addition  to  air  temperature,  the  air  pressure  has  an  effect  on 
I.D..  Figure  15  for  DF-2  shows  that  increasing  the  mean  pressure 
(during  I.D.)  reduces  the  I.D.  of  the  fuel  at  different  air  tempera¬ 
tures.  This  is  believed  to  be  mainly  due  to  the  effect  of  air  pressure 
on  the  preignition  chemical  reaction  rates.  This  effect  has  been 
reported  recently  by  Cowell  and  Lefebure  [26]  for  homogeneous  hydrocar¬ 
bon-air  mixtures  in  a  continuous  flow  apparatus.  The  gaseous  fuels  used 
are  propane,  ethylene,  methane  and  acetylene.  The  temperatures  ranged 
between  670°K  to  1020*K,  the  pressures  ranged  between  1  to  10  atmos- 
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pheres,  the  equivalence  ratios  ranged  between  0.2  to  0.7  and  velocities 
ranged  between  55  and  30  m/sec.  The  dependence  of  I.D.  on  pressure  has 
the  following  form 


I.D.  a  P 


where  n  varies  from  1  for  methane  to  0.66  for  acetylene.  For  liquid 
fuel  sprays  injected  in  diesel  engines  the  results  of  figure  15  show 
that  n  is  1.2  for  DF-2  (27). 

The  value  of  E  in  equation  (3)  can  be  obtained  by  plotting  log  I.D. 
vs  at  constant  Pm,  as  shown  in  figure  16  for  the  data  of  figure  15. 
The  slope  of  the  straight  line  in  figure  16  indicate  that  the  activation 
energy  E  for  DF-2,  at  Pm  =  400  psi  is  4965  BTU/16  mole.  Birdi  [27] 
found  that  E  is  a  function  of  the  mean  pressure  during  I.D.  in  addition 
to  the  C.N.  of  the  fuel.  Figure  17  shows  the  variation  in  E  with  the 
C.N.  at  Pm  =  400  psi,  for  all  the  reference  and  distillate  fuels  and 
their  blends,  tested  in  this  program,  at  ambient  air  temperatures  rang¬ 
ing  between  -4°F  (-20°C)  and  350°F  (176. 7°C). 

It  is  recommended  that  such  relationships  between  E  and  C.N.  be 
obtained  at  different  Pm.  Then  equation  (3)  would  be  able  to  correlate 
I.D.  with  C.N.  of  any  fuel. 

2.3.5  Inaccuracy  or  C.N.  in  Ranking  the  Overall  Performance  and 
Emission  Behaviour  of  Fuels  in  Diesel  Engines 

Many  recent  attempts  have  been  made  to  characterize  the  overall 
performance  and  emission  behaviour  of  fuels  in  diesel  engines  based  on 
C.N.  (21,22,28,29,30).  The  following  nine  performance  and  emissions 


32 


LOGARITHM  OF  IGNITION  DELAY  PERIOD  (MSEC) 


o 

DF2 

(CN-45.5) 

At 

RF1 

(CN-27.8X, 

X 

RF2 

(CN-30.0) 

❖ 

JP4 

(CN-35.0) 

z 

CN  - 

100 

V 

CN  - 

67 

X 

CN  - 

35  / 

* 

CN  - 

25  / 

z 

CN  - 

15  / 

* 

CN  - 

8  r 

1/TM  • 1 0**4  (1/R) 


Fig.  (16): 


Log  I.D.  plotted  vs  ^  for  fuels  and 
blends  of  different  C.N.S. 


100  FSXA-  0  DF2 ,  CN  2,  a 
[2.76  c  Rpi,  CN  =  27  ' 

«  RF2 ,  CN  -  1° 

O  JF“>  CN  *  3 
*  CM  =  1°° 

A  CN  *  67 

4.  CN  *  35 

n  cn  =  25 

«  CN  =  ^ 

ACH-  8 
v  CN  s  0 


ref 


CETAN& 


(17  V* 


jrisVilP 


\>eW< 


BeiaPi®  400  PSI* • 
at 


parameter  have  been  considered:  1)  Brake  thermal  efficiency,  2)  Noise, 
3)  Unburned  hydrocarbon  emissions,  4)  Nitrogen  oxide  emissions,  S) 


Carbon  monoxide  emission,  6)  Black  smoke  intensity,  7)  Particulate 
emissions,  8)  Startability,  and  9)  White  smoke. 

In  the  author's  opinion  it  is  impossible  to  develop  a  single  com¬ 
bustion  parameter  which  would  be  able  to  rate  the  performance  and  emis¬ 
sions  behaviour  of  fuels  in  diesel  engines  of  different  designs.  The 
reasons  behind  this  opinion  are: 

1.  Any  rating  requires  a  scale  with  arbitrary  reference  states 
representing  the  lowest  and  the  highest  ratings  on  the  scale. 
In  order  to  develop  a  scale  which  would  be  able  to  rate  the 
nine  performance  and  emission  parameters  two  fuels  are  needed: 
(i)  a  fuel  which  will  produce  the  highest  brake  thermal  effi¬ 
ciency,  lowest  noise,  lowest  HC,  CO,  N0X,  smoke,  particulates, 
and  white  smoke.  This  in  addition  to  be  the  easiest  to  start 
the  engine  on;  (ii)  similarly  another  fuel  would  be  needed 
which  would  produce  the  worst  performance.  It  is  clear  that 
it  is  impossible  to  find  these  two  reference  fuels  which  has 
the  extremes  in  all  the  nine  categories.  Probably  n-cetane 
does  not  represent  the  best  fuel  in  fuel/engine  performance 
parameters  nor  HMN  represents  the  worst  fuel. 

2.  Even  if  the  above  two  references  fuels  are  identified  it  is 
almost  impossible  that  different  fractions  in  their  blends 
would  give  corresponding  rating  on  the  scale  for  the  nine 
parameters.  Each  of  these  parameters  is  dependent  in  a  very 
complex  manner  on  many  independent  fuel,  air,  and  engine 
variables.  Even  most  of  the  relationships  between  many  per- 
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formance  parameters  and  the  independent  variables  are  not 
known  or  well  identified. 

The  C.N.  is  one  of  the  parameter  which  affect  the  I.D.  of  fuels  as 
explained  earlier.  Any  performance  parameter  which  is  a  strong  function 
of  I.D.  might  be  roughly  rated  by  the  C.N.  Doyle  and  Needhan  [20,21] 
reported  that  I.D.  defines  cold  start,  noise,  light  load  hydrocarbons, 
and  light  load  CO.  In  this  respect  C.N.  might  be  a  rough  indicator  of 
the  behaviour  of  fuels  in  engines.  Doyle  and  Needham  [20,21]  found  that 
the  peak  rate  of  heat  release  (PRHR)  correlated  well  with  NO^  emissions. 
PRHR  is  not  only  a  function  of  I.D.,  but  also  of  the  rate  of  mixing 
during  the  ignition  delay  and  the  premixed  combustion  periods ,  and  the 
rate  of  the  oxidation  reactions  during  the  premixed  combustion  period. 
Accordingly  C.N.  should  not  be  able  to  rate  NO^  at  all  operating  condi¬ 
tions  which  affect  the  premixed  combustion  rates. 

In  the  authors  opinion  an  approach  based  on  the  work  of  Wade  and 
Hunter  [30]  may  be  useful  in  correlating  the  fuel/engine  performance 
parameters  and  the  C.N..  Wade  and  Hunter  [30]  introduced  the  following 
three  combustion  parameters,  as  illustrated  in  figure  18: 

i.  Premixed  Combustion  Fraction  (PCF). 

ii.  Premixed  Combustion  Index  (PCI), 

iii.  Diffusion  Combustion  Index  (DCI). 

The  premixed  combustion  fraction  P,  is  the  fraction  of  the  total  heat 
released  at  the  end  of  the  premixed  combustion  mode.  The  premixed 
combustion  index  P/P^  is  the  fraction  of  the  total  energy  available  in 
the  fuel  present  in  the  cylinder  at  the  end  of  the  premixed  mode.  The 

diffusion  combustion  index  is  the  ratio  of  the  minimum  period  of  time 

9 

for  diffusion  combustion  Djj  to  the  actual  period  of  time  for  the  diffu- 
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sion  combustion,  0^.  A  diffusion  combustion  fraction  D  may  be  defined 
as  the  fraction  of  the  total  energy  released  in  the  diffusion  combustion 
mode.  Wade  and  Hunter  applied  their  proposed  technique  to  several  types 
of  diesel  engines  including  D.I.  and  I.D.I.  engines.  Two  fuels  having 
C.N.  =  47.5  and  C.N.  =  35  were  used.  The  effect  of  reducing  the  C.N.  on 
the  proposed  parameters  were  as  follows: 

i.  1.0.  increased. 

ii.  Both  the  premixed  combustion  fraction  and  the  premixed 
combustion  index  increased  by  more  than  100%  due  to  the 
long  I.D.  of  the  low  C.N.  fuel. 

iii.  The  diffusion  combustion  index  decreased  indicating  that 
the  burning  rate  of  the  low  C.N.  fuel  had  decreased. 

The  PCF  affects  HC,  and  N0^  and  combustion  noise,  and  the  DCI 
affects  BSFC  and  smoke. 

2*3*6  Problem  No.  6:  Inability  of  C.N.  in  Rating  the  Cold  Starta- 

bility  of  Diesel  Engines  on  Different  Fuels: 

The  difficulty  in  starting  diesel  engines  in  low  ambient  tempera¬ 
tures  represents  one  of  the  serious  concerns  of  diesel  engine  designers, 
particularly  when  alternate  fuels  are  used.  The  engine  startabil  ity 
depends  on  a  combination  of  factors  among  which  the  combustion  failure 
is  the  most  common.  Combustion  failure  during  starting  may  result  from 
improper  preignition  or  post  ignition  processes,  assuming  the  the  fuel 
flow  is  not  interrupted  at  the  low  ambient  temperatures.  The  preigni¬ 
tion  processes  control  the  length  of  I.D.  and  the  post  ignition  combus¬ 
tion  processes  control  the  rate  and  extent  of  the  oxidation  reactions. 
These  processes  are  affected  by  the  conditions  during  engine  starting, 
which  are  characterized  by: 
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'  i.  Low  ambient  air  temperatures  which  result  in  low  tempera- 

■  tures  near  the  end  of  the  compression  stroke. 

ii.  Low  engine  rotating  speeds  which  result  in 

(a)  low  compression  pressures  and  temperatures, 
•  (b)  poor  injection  process  resulting  in  poor  atomi¬ 

zation, 

^  (c)  low  air  turbulence  resulting  in  low  heat  and 

mass  transfer  rates  between  the  liquid  and  gas- 
^  eous  phases,  and  poor  mixing. 

The  fuel  properties  play  a  big  role  in  achieving  successful  igni¬ 
tion  and  subsequent  combustion.  The  most  important  parameter  to  achieve 
successful  ignition  is  the  length  of  the  ignition  delay  period.  High 
C.N.  fuels  have  short  ignition  delays  and  can  produce  successful  igni¬ 
tion.  Needham  and  Doyle  [21,22]  indicated  that  for  the  three  middle 

1 

distillate  fuels  tested  [1,2,  and  10]  cold  startability  correlates  well 
with  the  cetane  rating. 

Ignition  may  occur,  however  the  engine  may  not  start  because  of  the 
low  rate  of  post  ignition  oxidation  processes.  This  may  result  from  the 
lack  of  fuel  evaporation  or  mixing  the  vapor  with  the  air.  It  seems, 
then,  that  the  fuel  volatility  should  play  a  big  role  in  the  post  igni¬ 
tion  processes. 

Combustion  failure  tests  have  been  conducted  on  a  single  cylinder 
research  diesel  engines,  in  which  the  ambien*  temperature  is  lowered 
till  combustion  failed.  The  test  are  conducted  on  DF-2  (C.N.  =  45),  a 
light  fuel  (C.N.  =  27.8),  and  a  heavy  fuel  (C.N.  =  30).  The  properties 
these  fuels  are  given  in  Appendix  D. 
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Figure  19,20  and  21  show  the  cylinder  pressure,  needle  lift  and 
crankangle  degree  marks  for  different  ambient  temperatues.  The  traces 
are  for  0°C  and  lower  ambient  temperatures.  The  traces  for  higher 
ambient  temperatures,  in  the  range  between  0°C  and  the  normal  room  tem¬ 
perature,  are  not  included  in  the  figures  because  they  have  been  found 
to  be  very  close  in  shape  to  those  at  0°C  for  the  three  fuels  tested. 
Figure  19  is  for  DF-2,  figure  20  is  for  Reference  2  (Heaviest)  fuel,  and 
figure  21  is  for  Reference  1  (most  volatile)  fuel.  The  fuel  properties 
are  given  in  Appendix  D.  Figure  19  shows  that  all  the  traces  are  repro¬ 
ducible.  The  I.D.,  rate  of  pressure  rise  and  maximum  pressure  increase 
gradually  with  the  drop  in  ambient  temperature.  This  indicates  an 
increase  in  the  premixed  combustion  fraction. 

Figure  2®  is  for  Reference  2  fuel  and  shows  the  following: 

i.  I.D.  increases  with  the  drop  in  ambient  temperature, 
indicating  a  decrease  in  the  rates  of  the  preignition 
processes.  It  is  not  clear  at  this  time  whether  the 
physical  or  the  chemical  preignition  processes  has  con¬ 
tributed  to  the  increase  in  the  I.D.  because  the  fuel  is 
very  heavy. 

ii.  The  injection  process  appeared  not  to  suffer  as  the 
ambient  temperature  decreased.  The  fuel  pressure  and 
needle  lift  showed  no  cycle-to-cycle  variations  except  at 
an  ambient  temperatue  of  -25°C,  when  misfiring  took 
place. 

iii.  The  slope  of  the  pressure  rise  due  to  premixed  combustion 
decreased  gradually  with  the  drop  in  ambient  temperature, 
inspite  of  the  increase  in  I.D.  This  indicates  a  reduc¬ 
tion  in  premixed  combustion  index. 
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iv.  The  total  pressure  rise  due  to  combustion  decreased  with 
the  drop  in  temperatue.  This  indicates  a  drop  in  the 
total  heat  release. 

v.  Misfiring  is  being  caused  by  the  failure  of  the  premixed 
combustion  regime  which  resulted  in  the  failure  of  diffu¬ 
sion  combustion. 

Figure  21  is  for  Reference  1  (most  volatile  fuel)  and  shows  the 
following: 

i.  The  injection  process,  particularly  the  needle  lift 
showed  cycle-to-cycle  variations  at  a  temperature  of 
-20°C,  and  increased  with  further  drop  in  temperature, 

ii.  1.0.  increased  with  the  drop  in  ambient  temperatures, 

iii.  Cycle-to-cycle  variations  occurred  at  a  temperature  of 
-20 °C,  and  increased  with  further  drop  in  temperature. 
The  cause  of  this  variation  is  not  clear,  whether  it  is 
from  the  injection  or  the  preignition  chemical  processes. 
For  this  volatile  fuel  it  is  clear  that  the  preignition 
physical  processes  do  not  suffer  from  the  drop  in  tem¬ 
perature. 

iv.  The  rate  of  pressure  rise  due  to  premixed  combustion 
showed  cycle-to-cycle  variation,  corresponding  to  that  of 
the  I.D. 

v.  The  failure  of  combustion  cannot,  at  this  time,  be  attri¬ 
buted  to  the  failure  in  the  premixed  combustion  regime. 

The  above  analysis  of  the  experimental  results  on  combustion  fail¬ 
ure  at  900  RPM  for  the  three  fuels  indicates  the  following: 
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i.  The  volatility  has  a  great  effect  on  the  behaviour  of 
fuels  at  low  ambient  temperatures, 

ii.  The  failure  of  combustion  at  low  ambient  temperature  is 
caused  by  the  failure  of  the  premixed  combustion  for  the 
heavy  fuel. 

iii.  The  cause  of  the  failure  of  combustion  of  the  volatile 
fuel  is  not  very  clear,  because  of  the  deterioration  of 
the  injection  process.  However,  it  appears  that  the 
premixed  combustion  did  not  suffer  as  it  did  with  the 
heavier  fuel. 

Needham  and  Doyle  [21,22]  reported  similar  results  with  some  volatile 
fuels.  They  found  that  the  broadcut  fuel  (no.  4),  which  has  the  highest 
volatility  (distillation  range  64°C  to  312°C),  lowest  molecular  weight 
and  lowest  viscosity  of  all  the  fuels,  produced  the  worst  performance  in 
the  D.I.  engine  by  a  large  margin,  in  spite  of  a  similar  C.N.  to  that  of 
tar  sands  fuel  (no.  10). 

It  is  clear  that  the  cold  startability ,  like  the  other  fuel/engine 
parameters  depends  on  other  parameters,  in  addition  to  the  C.N.  of  the 
fuel . 

Montemayor  and  Owens  (1985)  determined  the  "Minimum  Unaided  Start¬ 
ing  Temperature"  for  fuels  of  different  physical  properties  and  cetane 
numbers.  They  reported  that  the  Minimum  Unaided  Starting  Temperature 
for  a  4  cylinder,  2-cycle  D.I.  diesel  engine  to  be  related  to  viscosity, 
50%  Boiling  Point,  and  an  autoignition  temperatue,  in  addition  to  C.N. 

The  problems  of  cold  startability  and  white  smoke  have  not  been 
well  researched  and  deserve  a  careful  investigation. 
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2.4  CONCLUSIONS 


1.  The  basic  function  of  the  Cetane  Scale  is  to  rate  the  autoignition 
behaviour  of  a  fuel  by  matching  its  ignition  delay  to  that  of  a 
blend  of  two  primary  Reference  fuels  under  specified  conditions  of 
inlet  temperature,  engine  speed,  coolant  temperature,  fuel-air 
ratio,  injection  timing  and  injector  opening  pressure.  The  only 
engine  variable  is  the  compression  ratio.  The  Cetane  Scale  was  not 
intended  to  rate  any  other  combustion,  performance  or  emissions 
parameter. 

2.  Three  intrinsic  problems  exist  in  the  cetane  rating  of  fuels. 

2.1  Inability  to  rate  low  ignition  quality  fuels.  The  solution  is 
to  raise  the  inlet  air  temperature  above  the  150°F  specified 
in  the  ASTM  D-613  procedure. 

2.2  Instability  of  the  ignitioa  delay  meter  in  rating  low  ignition 
quality  fuels.  This  has  been  found  to  be  due  to  the  cycle-to- 
cycle  variation  in  the  I.D.  of  these  fuels.  No  solution  is 
known  at  this  time. 

2.3  Inconformity  of  the  CFR  Combustion  System  with  Actual  Diesel 
Combustion  Systems.  A  possible  solution  is  to  redesign  the 
combustion  system.  Part  of  the  problem  may  be  solved  for  the 
existing  system  by  optimizing  the  injection  pressure  such  that 
the  rating  of  the  fuel  is  independent  of  the  needle  opening 
pressure . 

3.  Three  fuel/engine  preformance  problems  are  examined: 

3.1  Invalidity  of  C.N.  in  correlating  the  I.D.  for  different  fuels 
in  different  engines  under  all  operating  conditions.  The  I.D. 
should  be  correlated  in  terms  of  three  basic  parameters:  mean 
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temperature,  mean  pressure,  and  a  global  activation  energy. 
The  global  activation  energy  can  be  correlated  with  the  C.N. 


3.2  Inaccuracy  of  C.N.  in  rating  the  overall  performance  and 

\ 

emissions  behaviour  of  fuels  in  diesel  engines.  This  correla¬ 
tion  should  be  done  in  terms  of  I.D.  and  the  combustion  para¬ 
meters  which  characterize  the  premixed  and  diffusion  combus- 
>  tion  regimes. 

3.3  Inability  of  C.N.  in  rating  the  cold  startability  of  diesel 
engines  on  different  fuels.  Startability  depends  on  the 
success  of  both  the  autoignition  and  the  subsquent  combustion 
reactions.  Many  parameters  other  than  C.N.  need  to  be 
included. 

2.5  RECOMMENDATIONS 

1.  Study  the  cetane  rating  at  inlet  air  temperatures  higher  than  the 
150°F  specified  by  the  ASTM  D-613  procedure.  Determine  the  most 
suitable  temperature  which  allows  a  suitable  size  of  electric 
heater  to  be  used. 

2.  Invistigate  the  phenomena  of  cycle-to-cycle  variation  in  the  auto¬ 
ignition  process  of  fuels. 

3.  Investigate  the  effect  of  optimizing  the  needle  opening  pressure 
for  each  fuel  on  the  cetane  rating. 

4.  Develop  correlations  between  the  C.N.-  and  the  global  activation 
energy  for  the  preignition  reactions  for  different  fuels,  to  be 
used  for  correlating  I.D.  with  C.N.,  temperatue  and  pressure. 

5.  Develop  correlations  between  I.D.  and  the  combustion  parameters  for 
the  premixed  and  diffusion  combustion  and  the  fuel/engine  perfor- 
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mance  parameters.  Relate  these  parameters  to  C.N.  by  using  C.N.- 
I.D.  correlations. 

6.  Make  a  combustion  failure  analysis  during  engine  starting  on  diffe¬ 
rent  fuels  and  ambient  temperatures.  This  will  identify  the  engine 
and  fuel  parameters,  other  than  C.N.,  which  affect  engine  starta- 
bility. 
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2.6  APPENDICES 


2.6. 1  APPENDIX  A 


CFR  ENGINE  SPECIFICATIONS 


Bore 

Stroke 

Displacement 
Compression  Ratio 
Connecting-rod,  center-to-center 
Number  of  piston  rings 
Crankcase 
Injection  Pump 


Injector 

Injector  Pintle  valve  lift 
Valve  timing 

Intake  opens 
Intake  closes 
Exhaust  opens 
Exhaust  closes 


3.25  in. 

4.50  in. 

37.33  cu  in. 

7  to  28 
10.00  in. 

5 

Highspeed  CFR 
Bosch  Type  of  vari¬ 
able  injection 
timing 

Bosch,  ADN-30S-3/1 
.004  to  .006  in. 

10°  ATDC 
34°  ABDC 
40°  BBDC 
15°  ATDC 
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2.6.2  APPENDIX  B 


SUMMARY  OF  OPERATING  CONDITIONS  SPECIFIED 
BY  ASTM-D-613  CETANE  METHOD 


Engine  speed,  rpm 

Crankcase  oil,  SAE  grade 

Oil  pressure  at  operating  temperature 

Crankcase  oil  temperature 

Coolant  Temperature: 

Range 

Constant  within 
Intake  Air  Temperature 
Valve  clearances 


900  ±  9 
30 

25  to  30  psi 
135  ±  15°F 

212  ±  3°F 
±  1°F 
150  ±  1°F 

.008  in.,  measured 
with  the  engine  hot 
and  running  under 
standard  operating 
conditions  on  a 
reference  fuel  of  50 


C.N. 


Fuel  Injection  Pump,  Bosch: 
Advance,  deg  btdc 
Opening  pressure,  psi 
Injector: 

Water  jacket  temperature 
Fuel  injection  rate 


13.0 

1500  ±  50 

100  t  5°F 

13.0  t  .2  ml/min. 


2.6.3  APPENDIX  C 
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VISCOSITY  (cSt  AQ°C)  2. AO  2.76  A. 37  0.91  33.93 


Properties  of  DOE  Fuels  Tested  by  Reference  [21.22 

(continued) 

(NOS.  6  through  10) 
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HEAT  OF  COMBUSTION 

(NET  MJ/kg)  42.61  42.21  42.25  43.36  37.03 


Properties  of  DOE  Fuels  Tested  by  Reference  [21,22} 
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HEAT  OF  COMBUSTION 

(NET  MJ/kg)  37.23  43.07  41.21  39.72  41.98 


Properties  of  Fuels  Tested  at  WSU 
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SECTION  3 


C.N.  RATING  OF  COLD  STARTABILITY  OF  DIESEL  ENGINES 
ON  DIFFERENT  FUELS 


The  inability  of  C.N.  scale  in  rating  the  cold  startability 
of  diesel  engines  on  different  fuels  is  of  great  concern  to  the 
military  and  industry.  This  seciton  of  the  report  deals  with 
the  research  conducted  to  find  out  why  the  C.N.  fails  to  predict 
the  cold  startability  of  diesel  engines. 

This  section  describes  the  experimental  work  conducted  on 
a  single  cylinder  engine  to  determine  the  effect  of  the  Cetane 
number  and  other  fuel  properties  on  cold  startability. 

3.1  ENGINE  AND  INSTRUMENT ION : 

The  engine  is  a  single  cylinder,  Deutz,  4-stroke-cycle, 
direct  injection  air  cooled  diesel  engine  d  =  0.095  m,  and  C.R. 

=  17:1. 

The  cylinder  gas  pressure  is  measured  by  using  a  quartz 
crystal  pressure  transducer,  and  a  charge  amplifier.  Before 
eacu  run  the  engine  was  soaked  at  the  controlled  ambient 
temperature  for  at  least  5  hours. 

The  instantaneous  angular  velocity  Is  determined  from  the 
time  taken  by  the  crankshaft  to  rotate  through  a  fixed  angle. 

An  electromagnetic  pick  up  installed  in  the  flywheel  housing, 
produces  a  pulse  signal  to  indicate  the  present  of  a  tooth  on 
the  flywheel  ring  gear. 

The  data  acquisition  system  consists  of  a  24  bit  high  speed 
counter,  four  10  bit  analog  to  digital  converters  and  a  TDC 
Identification  circuit. 

3.2  FUELS: 

The  following  two  fuels  are  used:  DF-2  (Diesel  fuel  r umber 
2)  and  JP-4  fuel.  The  fuel  main  specifications  are  given  in 


item  2.6.4*  JP-4  is  a  lighter  fuel  than  DF-2,  and  has  a  lower 
C.N.,  but  a  higher  volatility  than  DF-2. 

3.3  EXPERIMENTAL  RESULTS 

A.  DF-2  FUEL 

Fig. 22  is  a  plot  of  the  instantaneous  angular  velocity  of 
the  engine  versus  the  engine  revolutions,  at  an  ambient  temperature 
of  20°C,  after  a  cold  soak  period  of  five  hours.  The  engine  started 
successfully  and  reached  the  idling  speed  of  130  radians  per  second 
after  about  30  revolutions.  Fig  23  is  for  starting  the  engine  at 
an  ambient  Temperature  of  15°C,  after  a  cold  soak  period  of  fire 
hours.  The  data  for  20°C  is  plotted  on  Fig  22at  the  same  scale  for 
comparison.  The  engine  started  successfully  and  reached  the  idling 
speed  after  35  revolutions.  The  oscillations  in  the  engine  speed 
after  the  engine  starts  are  under  investigation  at  the  present  time. 
Fig  24  is  at  an  ambient  temperature  of  10°C  after  a  cold  soak  time 
12  hours.  The  engine  started  and  reach  the  idling  speed  after  45 
revolutions.  The  oscillations  in  speed  are  more  frequent  and  their 
amplitude  is  higher  than  at  the  higher  ambient  temperatures.  Fig  25 
is  for  an  ambient  temperature  of  5°C  and  cold  soak  of  5.5  hours. 

The  uniform  cranking  speed  is  60  radians  per  second,  for  about  45 
revolutions.  Following  this  the  engine  began  to  accelerated  at  a 
very  slow  rate,  which  increased  in  time.  A  maximum  speed  of  138 
radians  per  second  is  reached  after  90  revolutions  from  the  start 
of  cranking,  or  after  45  revolutions  from  uniform  cranking.  Fig  25 
shows  that  the  oscillations  increased  in  number  and  amplitude  as  the 
temperature  decreased  from  10°C  to  5°C.  Fig  26  is  for  an  ambient 
temperature  of  2°C.  The  engine  is  cranked  at  a  uniform  speed  of  57 
radians  per  second,  for  about  60  revolutions  after  which  the  speed 
increased  very  gradually,  at  a  slow  rate  at  fiist.  The  rate  of 
increase  in  speed  increased  till  it  reached  a  maximum  of  143  radians 
per  second.  The  number  of  revolutions  taken  for  the  engine  to  accele¬ 
rate  from  the  uniform  cranking  speed  to  the  maximum  speed  is  50. 

The  oscillations  in  speed  increased  in  number,  and  amplitude. 
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REVOLUTIONS  OF  THE  ENGINE 


Plot  of  angular  velocity  vs.  engine  revolutions.  Fuel 
DF  2,  Oil:  arctic  oil,  ambient  temperature : 20  C,  cold 
soak  time:  5  hrs . ,  startability  conditions  Started- 
easily . 
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easily . 
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30  Plot  of  angular  velocity  vs.  crank  angle  ,  fuel:  JP  4 

Bilrarctlc  oil,  cold  soak  time:  12  hr»a,  ambient  tempe¬ 
rature:  25  C,  startability  condition:  started  easily. 


Fig.  31 


Plot  of  angular  velocity  vs.  crank  angle,  fuel:  JP4, 
oil:  arctic  oil,  cold  soak  time:  6  hrs.,  ambient 
temperature : 10  C,  startability  condition:  started  easily 
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Fig  27  is  for  an  ambient  temperature  of  -5°C  and  Fig28is  for  an 
ambient  temperaute  of  -S°C.  By  comparing  the  two  figure  it  is  clear 
that  the  drop  in  ambient  temperature  resulted  in  a  decrease  in 
the  cranking  speed,  an  increase  in  the  cranking  period,  and  an  in¬ 
crease  in  the  amplitude  of  the  oscillations.  Fig  29  is  for  an 
ambient  temperature  of  -10°C.  The  engine  failed  to  start  even  after 
cranking  for  675  revolutions.  Fig  30 is  for  engine  startability  with 
JP-4  at  an  ambient  of  25°C  and  cold  soak  time  of  12  hours.  The 
engine  started  quickly  and  accelerated  to  137  radians  per  second  in 
20  revolutions.  Fig  31  is  for  JP-4  and  an  ambient  temperature  of 
10°C.  Fig  32  is  for  an  ambient  temperature  of  5°C.  Fig  33  is  for 
an  ambient  temperature  of  0°C.  Fig  34  is  for  an  ambient  temperature 
of  -5°C.  Fig  35  is  for  an  ambient  temperature  of  -6°C.  Fig  36  is 
for  an  ambient  temperature  of  -8°C.  Fig  37  is  for  an  ambient  tem¬ 
perature  of  -9°C.  Fig  38  is  for  an  ambient  temperature  of  -10°C. 

Fig  39  is  for  an  ambient  temperature  of  -11°C.  Fig  40  is  for  an 
ambient  temperature  of  -13°C 

COMPARISON  BETWEEN  STARTABILITY  ON 
DP-2  AND  JPA: 


Fig  41  is  a  plot  of  the  engine  revolutions  from  the  beginning  of 
cranking  to  the  time  the  engine  reached  its  maximum  starting  speed, 
at  different  ambient  temperatures.  The  difference  between  DF-2 
and  JP-4  is  negligible  at  normal  ambient  temperatures.  As  the 
ambient  temperature  drops  the  difference  between  the  two  fuels  in¬ 
crease?.  The  engine  failed  to  start  on  DF-2  after  cranking  the 
engine  for  675  revolutions,  at  an  ambient  temperature  of  -10°C. 

At  the  same  ambient  temperature  the  engine  successfully  started  on 
JP-4  after  cranking  for  270  revolutions.  The  engine  did  not  fail 
to  start  on  JP-4  at  ambient  temperatures  down  to  -13°C. 
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Fig.  41  Comparison  between  the  engine  revolution  to 

maximum  starting  speed  on  DF  2  and  JP4  fuels. 
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It  is  clear  that  JP-4  is  more  volatile  than  DF-2.  However  the 
Cetane  number  of  JP-4  is  less  than  DF-2. 

The  experimental  results  indicate  very  clearly  that  the  fuel 
volatility  plays  a  major  role  in  the  formation  of  the  combustible 
mixture  which  is  ready  for  autoignition  at  the  low  ambient  tempera¬ 
tures.  Since  the  Cetane  number  is  determined,  according  to  ASTM 
D-613  method,  at  hot  air  and  engine  temperatures  it  does  not  re¬ 
flect  the  effect  of  fuel  volatility. 

It  can  be  concluded  that  the  C.N.  alone  cannot  rate  the 
autoignition  characteristics  of  the  fuel  at  air  and  engine  temr 
peratures  lower  than  those  used  in  the  ASTM  D-613  method.  These 
temperatures  are  150°F  for  inlet  air  temperature  and  220°F  for  the 
engine  coolant  temperature.  These  temperatures  do  not  represent 
®t  all  the  conditions  during  cold  starting. 
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The  effect  of  fuel  volatility  on  the  ignition  delay 
was  investigated  in  a  single  cylinder  research  diesel 
engine,  under  steady  state  running  conditions.  The 
tests  were  conducted  at  low  room  temperatures.  The 
ignition  delay  was  measured.  The  ambient  temperature 
was  lowered  till  the  engine  misfired.  A  comparison  was 
made  between  the  misfiring  temperatures  and  ignition 
delays  of  fuels  having  equal  cetane  numbers  but  of 
different  volatilities.  The  low  volatility  fuels  were 
composed  by  blending  primary  reference  fuels  so  that 
they  have  cetane  numbers  equal  to  thoxeof  distillate 
fuels.  The  primary  reference  fuels  are  cetane  (C16H34) 
and  heptamethyl nonane  (C15H34) .  Cetane  has  a  C.N.  of 
100  and  heptamethylnonane  has  a  C.N.  of  15. 

4.1  ENGINE 


The  engine  is  the  ATAC .  Hi . LABE CO  single  cylinder, 
direct  injection  diesel  engine,  4. 5x4. 5  in,  C.R.=16:1. 
The  engine  is  coupled  to  a  General  Electric  D.C. 
*fotor/Generator  Cradle  type  dynamometer.  The  dynamome 
ter  is  provided  with  a  speed  and  load  control  and  a 
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tachogenerator  for  measuring  the  speed .  The  intake  and 
exhuast  sides  were  connected  to  two  steel  surge  tanks 
of  3  cubic  feet  capacity  each  in  order  to  reduce 
pulsations  on  the  measuring  system. 

4.2  INSTRUMENTATION 

Figure  **2  shows  the  layout  of  the  equipments  and 
general  instrumentations.  Intake  air  is  taken  from  the 
cold  room,  which  has  vents  to  allow  for  the  make-up  air 
to  flow  into  the  room  and  be  cooled  before  it  reaches 
the  intake  air  filter.  A  butterfly  valve,  remotely 
controlled,  is  used  to  control  the  air  pressure  in  the 
intake  system.  The  air  pressure  is  adjusted  so  that 
the  mean  pressure  during  the  ignition  delay  period  is 
kept  constant  for  all  the  runs.  The  flow  rate  was 
measured  by  a  Meriam  Laminar  Flowmeter  (50MW20-2) , 
which  is  installed  before  the  intake  surge  tank.  The 
pressure  differential  across  the  flowmeter  is  measured 
by  an  inclined  manometer  with  a  scale  division  of  0.01 
inch.  The  flowmeter  was  calibrated  by  the  manufac¬ 
turer.  A  Kistler  (601B-SN  2922)  quartz  pressure 
transducer  is  installed  in  the  cylinder  head  of  the 
engine  to  measure  the  cylinder  pressure.  The  signal 
from  the  transducer  is  amplified  by  a  Kistler  504E 
charge  amplifier.  This  signal  is  then  fed  to  the 
voltage  channel  of  the  Textronix  565  oscilloscope,  and 
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<gure  42  Layout  of  Experimental  Setup 


fed  to  a  peak  meter,  in  order  to  measure  the  peak 
cylinder  pressure.  The  pressure  transducer  was 
statically  calibrated  using  a  dead  weight  tester 
machine,  over  a  range  0-1500  psi. 

Fig  43  shows  the  schematic  of  the  electronic 
instrumentation  used,  to  measure  the  needle  lift  of  the 
injector,  the  triggering  circuit,  the  injection 
pressure  transducers  and  the  pickup  used  to  indicate 
the  crank  angle  degrees. 

4.3.1  MISFIRING  TEMPERATURE 

The  misfiring  temperatures  for  distillates  and 
blend  of  primary  reference  fuels  are: 

DF2  (CN=45.5,  distillate)  :  no  misfiring  was 

found  at  -32°c 


Ref.  1  (CN-27.8,  distillate)  :  -28°c 

Ref.  2  (CN=30.0,  distillate)  :  -26°c 

JP4  (CN-35.0,  distillate)  :  -27°c 

CN-45.5  (blend  )  :  -9°c 
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'9ure43  Schematic  of  Electronic  Instrumentation  (Reference  12) 
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Figure  44'  shows  the  above  results  and  indicates 
that  the  misfiring  temperature  of  distillates  is  lower 
than  that  of  the  fuel  blends.  Because  all  the  fuel 
blends  in  this  research  work  were  formed  from  N-Cetane 
(CN-loo)  and  Heptamethylnonane  (HMN,  CN-15) ,  they  are 
less  volatile  than  the  distillates.  This  difference  in 
volatility  causes  the  large  difference  in  their 
misfiring  temperatures. 

4.3.2  IGNITION  DELAY 

The  effect  of  inlet  air  temperature  on  the  ignition 
delay  period  of  distillates  and  blends  of  primary 
reference  fuels  was  studied.  Figure  45  shows  that  the 
ignition  delay  periods  of  blended  fuels  are  higher  than 
those  of  distillated  fuels,  while  they  have  the  same 
cetane  number. 

Figure  46  shows  variation  of  ignition  delay  with 
cetane  number  at  different  ambient  air  temperatures. 
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MISFIRING  TEMPERATURE 
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20  30  40  50 
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Figure  ^6  Variation  of  ignition  Oelay  Period  with  Cetane  Number  at 


Different  Ambient  Temperature  (d i at i II ated  and  blend  fuels) 
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The  ignition  delay  period  for  distillates  are  lower 
than  the  ignition  delay  for  the  blends,  both  having  the 
same  cetane  number. 

4.4  CONCLUSIONS 

This  analysis  clearly  indicates  that  the  cetane 
number  of  the  fuel  is  not  the  only  parameter  which 
affect  the  autoignition  process  in  diesel  engines, 
particularly  at  low  ambient  temperatures.  The  shorter 
delay  of  the  distillates  and  lower  their  misfiring 
temperature  compared  with  the  blends  of  primary 
reference  fuels,  both  having  the  cetane  number,  support 
the  above  conclusion. 
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SECTION  5 


MATHEMATICAL  MODEL  FOR  ENGINE  STARTABILITY 


A  mathematical  model  has  been  developed  to  predict  the  behavior  of 
the  engine  during  the  cranking  and  starting  process.  The  model  will  be 
presented  in  a  paper  in  the  SAE  Annual  Congress  to  be  held  in  Detroit,  Mi. 
February,  19^8.  The  paper  is  duplicated  in  this  section. 


DIESEL  STARTING: 
A  MATHEMATICAL  MODEL 


Timothy  P.  Gardner 

Ford  Motor  Company 

Naeim  A.  Henein 

Wayna  Stata  Univarsity 


ABSTRACT 

A  mathematical  nodal  la  developed  co  study 
the  transient  behavior  of  a  four-stroke,  single 
cylinder  naturally-aspirated,  DI  diesel  engine 
during  cranking  and  starting.  The  model  simulates 
che  full  the rao dynamic  cycle  of  the  engine  and 
Includes  detailed  sub -mode Is  for  che  intake  and 
exhaust  gas  flow  processes,  autoignition 
combustion,  heat  transfer,  mechanical  friction, 
blowby,  and  engine  dynamics.  The  model  considers 
the  period  of  time  from  starter  on  until  che 
engine  reaches  che  idle  speed. 

Experimental  data  and  results  obtained  from 
Che  engine  are  used  co  callbraca  and  validate  Che 
model.  A  comparison  of  results  from  Che  modal  and 
che  experiments  generally  show  good  agreement  for 
Che  starting  conditions  investigated. 


MART  FACTORS  CONTRIBUTE  Co  che  successful 
star'  j.ng  of  diesel  engines.  In  che  absence  of  a 
positive  Ignition  source  In  che  diesel  engine, 
che  initiation  of  combustion  is  mainly 
temperature  dependent.  The  air  temperature  near 
che  end  of  compression,  In  turn,  depends  on  the 
compression  ratio  and  other  engine  operating 
conditions  such  as  inlet  air  temperature  and 
cranking  speed.  During  cranking,  che  compressed 
air  temperatures  are  lower  chan  chose  during  any 
ocher  mode  of  engine  operation.  The  low  compres¬ 
sion  temperatures  are  caused  partly  by  the  low 
ambient  temperature,  and  more  importantly  by  che 
excessive  heat  losses  between  the  cylinder  gases 
and  che  walls  and  by  che  excessive  blowby  losses 
at  Che  low  cranking  speeds.  Consequently,  the  low 
compression  temperatures  result  in  poor  ignitl- 
blllty  of  the  fuel,  and  hence,  poor  startabllicy 
of  the  engine. 

Host  engine  designers  agree  chat  scarclng 
requirements  set  the  lower  limit  on  che  range  of 
compression  ratios  used  in  diesel 
engines.  Therefore,  a  model  which  could  predict 
che  minimum  CR  at  che  design  stage,  which  would 
eventually  lead  to  a  successful  start,  would  be 


of  considerable  interest  to  che  designer. 

This  paper  describes  a  general  machematlcal 
model  of  a  naturally- aspirated ,  direct-injection 
(DI)  diesel  engine  which  can  be  used  co 
investigate  the  effects  of  compression  ratio, 
ambient  temperature,  fuelling  rate,  etc.  on  che 
starting  process.  The  model  simulates  a  full 
thermodynamic  cycle  In  which  cylinder  pressure, 
temperature,  equivalence  ratio,  etc.  are 
calculated  along  with  engine  friction  co 
determine  the  final  compression  temperatures  In 
the  cylinder  as  well  as  the  engine  cranking  speed 
during  starting.  The  model  also  calculates  the 
instantaneous  heat  losses,  blowby  losses,  end 
frictional  losses  during  che  starting  process. 

LITERATURE  REVIEW 

Many  machematlcal  models  have  been  developed 
co  simulate  che  fluid  mechanics  and  thermodynamic 
processes  in  Internal  combustion  engines  [1-11]. 
These  models  can  be  grouped  into  two  broad 
categories;  1>  steady-state  models,  and  2) 
unsteady-state  or  transient  models. 

In  steady- state  models,  che  engine  speed, 
fuel  race,  and  load  are  usually  held  constant 
over  che  entire  cycle.  By  definition,  steady 
state  operation  requires  that  the  conditions  at 
the  end  of  the  cycle  must  equal  Chose  at  che 
beginning.  This  criteria  Is  used  co  obtain  an 
energy  balance. 

In  transient  models,  thermodynamic  end  gas 
flow  processes  are  usually  simulated  by  using  one 
of  two  approaches.  The  first  is  a  quasl-linear 
approach,  ,ln  which  empirical  relationships  are 
used  for  engine  thermodynamics  and  gas  flow. 
These  are  linked  with  dynamic  models  of  the 
mechanical  components  [10].  The  major 
disadvantage  of  quasi- linear  models  Is  their 
heavy  reliance  on  experimental  data  which  limits 
che  usefulness  of  this  approach.  The  second  Is  a 
quasi-steady  approach,  in  which  engine  thermo  and 
fluid  dynamics  are  treated  on  a  quasi-steady 
basis  employing  che  'filling  and  emptying* 
concept  [4,7],  By  using  the  quasl-stsady 
approach,  one  can  calculate  the  steady,  or 


transient  phenomena  securing  in  Che  engine 
regardless— of  its  operating  condi clone.  However, 
the  formulation  end  soludon  procedure  must  be 
different  Chen  chee  used  in  steady  models  since 
che  periodicity  condition  (constenc  speed)  does 
not  hold  during  transient  operation. 

Most  filling  and  emptying  models  for  diesel 
transient  analysis  reported  in  che  literature 
have  been  used  Co  calculate  Che  performance 
characteristics  of  turbocharged  diesel  engines 
under  changing  speed  or  load,  or  to  match  che 
Curbocharger  to  a  new  or  uprated  engine 
design.  The  authors  know  of  no  one  who  has 
applied  the  filling  and  emptying  models  to  study 
the  transient  behavior  of  diesel  engines  during 
initial  cranking  and  starting.  Since  che  basic 
assumptions  used  in  quasi-steady  models  are 
applicable  to  -  any  operating  condition  in  the 
engine,  an  attempt  is  made  to  use  these  models  to 
develop  a  general  model  for  engine  startebllity 
and  to  investigate  che  effects  of  various  design 
and  operating  parameters  on  the  starting  process 
of  diesel  engines. 

The  primary  difference  between  che  present 
filling  and  emptying  model  and  some  of  che  ocher 
models  published  in  the  literature  is  in  che 
decailed  treatment  of  che  friction  losses  and 
blowby  losses,  which  are  particularly  important 
during  che  cranking  and  starting  process. 
Instantaneous  values  of  friction  and  blowby 
losses,  rather  than  cycle  averaged  values,  are 
calculated  along  with  instantaneous  heat  losses. 
The  Inclusion  of  these  calculations  In  che 
filling  and  emptying  model  allows  a  more  decailed 
Investigation  of  che  effects  of  various  design 
and  operating  parameters  on  these  losses,  as  well 
as  provide  additional  insights  Into  how  these 
parameters  affect  engine  startability . 

CRITERIA  FOR  ENGINE  STARTABILITY 

The  fact  Chat  ignition  occurs  in  the 
combustion  chamber  does  not  necessarily  mean  chat 
che  engine  will  successfully  start.  Marginal 
conditions  in  the  combustion  chamber  may  result 
in  partial  ignition  or  even  misfire,  which  may 
prevent  che  engine  from  starting.  Therefore,  a 
criteria  must  be  established  in  Che  model  for 
evaluating  che  startability  of  che  engine. 

A  number  of  different  methods  have  been  used 
in  experimental  studies  to  evaluate  engine 
startability.  These  methods  are  based  on;  1) 
required  cranking  time,  2)  time  co  start  of 
combustion,  3)  cranking  torque  &  smoke 
measurements,  4)  tlme-to-scarcer  off,  and  3) 
cime-co- idle  speed.  The  authors  believe  that  che 
clme-to-ldle  speed  method  is  che  most  reliable, 
and  the  most  repeatable  of  these  methods,  and 
hence,  is  adopted  for  che  starting  criteria  in 
che  present  model. 

MODEL  FOR  STARTING 

The  starting  model  used  in  this  study  is 
shown  schematically  in  Figure  1.  Initially,  the 
engine  is  assumed  co  be  sc  rest.  Starting 
consists  of  cranking  the  engine  while  fuelling 


Figure  1  -  Model  for  Starting  Process 


with  che  scarcer  motor  on  until  igniclon  occurs 
and  the  engine  accelerates  to  an  arbitrarily  sec 
cranking  speed.  The  starter  motor  torque  is  then 
Curned  off,  and  the  engine  continues  co 
accelerate  with  the  aid  of  Che  gas  pressure 
torque  until  it  reaches  its  idling  speed  to 
complete  che  starting  process.  The  starting 
attempt  is  considered  to  be  successful  if 
ignition  occurs  while  cranking  with  fuelling,  and 
Che  engine  is  able  to  accelerate,  with  no 
assistance  from  che  starting  motor,  co  its  idle 
speed  within  a  specified  time.  If  ignition  does 
not  occur  and/or  the  engine  falls  to  reach  its 
idle  speed  within  a  specified  time  (or  number  of 
cycles),  the  starting  attempt  is  considered 
unsuccessful. 

MATHEMATICAL  FORMULATION 

The  mathematical  formulation  of  che  model  is 
derived  by  applying  Che  thermodynamic  principles 
to  Che  three  control  volumes  shown  in  Figure  2 ; 
1)  che  cylinder,  2)  che  intake  manifold,  and  3) 
che  exhaust  manifold.  Heat,  work  and  mass 
transfer  across  the  control  surfaces  of  these 
control  volumes  are  calculated.  The  basic 
approach  used  to  derive  che  governing  equations 
for  each  control  voluma  is  to  apply  an  energy 
balance,  chan  combine  che  thermodynamic  property 
relationships,  and  che  mass  continuity  equation 
to  obtain  thi  final  equations. 

The  rate  of  change  of  gas  temperature  wich 
respect  co  Cima  dT/dt,  can  be  expressed  in  che 
form  of  an  ordinary  differencial  aquation 
developed  by  Borman  [3]  as; 
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where. 
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V  dt 

dt 
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■  it  *  - 

dc  dt 

M 

T  3R 

r  an 

P  3R 

D  -  1  - - 

R  3P 

Equation  (l)  la  cha  baalc  fora  of  cha  energy 
aquaclon  that  la  applied  Co  each  control  volume 
ac  aach  lnatane,  together  with  cha  continuity 
aquaclon,  and  cha  perfect  gaa  aquaclon.  All  of 
cha  partial  derivatives  in  Eq.(l)  can  avaluacad 
froa  cha  gaa  proparty  relationships  which  are  of 
cha  fora: 

u  -  u<  P,  T,  ♦  ) 

R  -  R(  P,  T,  «  ) 


All  valuaa  of  enthalpy  and  Internal  y 

used  In  the  nodal  ara  defined  on  cha  basis  of  c..,* 
absolute  reference  state  using  cha  definition  for 
absolute  enthalpies  given  by  Powell  [34].  The 
absolute  enthalpy  of  the  pure  fuel  Is  calcu)Med 
froa  the  known  enthalpy  of  coabustlon  and  Die 
sensible  enthalpy  of  the  fuel  according  to: 

hf  -  hf  ♦  Cp  dt  (3) 

where , 

hf  -  Che  enthalpy  of  foraation  froa  che 
aleaants  (at  T-0  °K  .and  P-1  a cm) 

Cp  -  che  specific  heat  of  che  fuel. 

The  voluaa,  V,  and  rate-of -change  of  volume, 
dV/dt,  for  the  cylinder  can  be  calculated  froa 
che  engine  geoaetry  and  speed  according  to  [37]  : 

V  -  vc  +  *D2S/4  (4) 

S  -  r  [1  -  cos®  +  L/r  (1-  A  -  r  sin2«/L2  ] 


The  rata-of-chango  of  aass  in  cha  cylinder, 
dM/dc,  obtained  by  considering  cha  races  of  flow 
through  che  Intake  and  exhaust  valve,  fuel 
Injection,  and  blowby  past  che  piston  rings: 


dV  v  N  D2  dS 


dt  120  d® 

N  D2 


1  cos® 


dM 

J*1 

«e 

—  + 

dMf 

dMb 

(2) 

where. 

-  r  sin®  l  +  - - - —  (5) 

120  (l  l2sin2«)  5 

dt 

dt 

dt 

dt 

dt 
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vc 

D 

-  the  clearance  voluaa  at  TOC 

-  che  cylinder  bore 

• 

A 


* 


w 


X  -  che  radon  r/L 
r  -  the  crank  throw  radius 
L  -  eha  connacdng  rod  length 
9  -  tha  crank  angle  measured  from  TDC. 

Tha  rata  of  haat  transfer,  dQ/dt,  at  tha 
gas -to -wall  Interface  can  ba  calculated  by: 

dQ 

-  -  hc  At  (  T-  -  Tw)  (6) 

dt 


cylinder  and  the  intaka  and  exhaust  aanlfolds. 
Tha  manifolds  are  treated  as  infinite  plenums 
with  specified  pressures,  thus,  neglecting  any 
dynamic  effects  in  the  flow  due  to  pressure 
waves.  From  a  knowledge  of  Che  gas  conditions  in 
the  cylinder  and  manifolds,  plus  a  geometric 
valve  area  and  discharge  coefficient,  che 
instantaneous  mass  flow  through  che  valves  (and 
ports)  can  be  calculated  by: 

«  *  Cd  Ae  Pt  (  Z  gc/R!  TjJ0-5  (11) 


where, 

he  -  che  gas -to -wall  heat  transfer 
coefficient 

At  -  che  total  wall  area  for  heat  transfer 

Tg  -  che  gas  temperature 

-  the  cylinder  wall  temperature . 

Calculation  of  che  heat  transfer  coefficient 
(hc)  for  che  cylinder  is  discussed  later. 

The  equivalence  ratio  •  in  the  cylinder  will 
change  whenever  mass  enters  or  leaves  che 
cylinder.  The  equation  for  9  at  any  time  t  can 
be  written  in  cents  of  the  mass  of  air  (Ma) ,  mass 
of  fuel  (Mf) ,  and  che  stociometrlc  fuel/air  ratio 
(f,)  as: 


Mf 

«  -  -  (7) 

«e  £| 

where, 

M 

Ma  -  -  (8) 

l  ♦  9fs 


M  9f„ 

Mf  -  -  (9) 

l  *  9f# 

Substitution  of  Eqs.(8)  and  (9)  into  Eq.(7),  and 
differentiating  Che  results  gives: 


d9  Fi 
dt  M 
where , 

Fl  *  1  ♦  ♦*. 

By  using  Eqs.(l).  (2),  and  (10)  to  evaluate 
dT/dt,  dM/dt,  and  d*/dc ,  new  values  of  T,  M,  and 
9  can  be  calculated  ac  Che  end  of  each  crank 
angle  step.  The  pressure,  P,  can  then  be 
calculated  using  Che  perfect  gas  equation. 

THERHOOYMAHIC  GAS  MXZOTRE  PROPERTIES  -  The 
equilibrium  thermodynamic  properties  for  the 
products  of  combustion  of  air  and  CnH2„  are 
assumed.  Algebraic  expressions  derived  by  Krieger 
and  Borman  [13]  are  used  for  lean  mixtures,  and 
expressions  by  Marxouk  [ 35 1  for  rich  mixtures. 

GAS  PROCESSES  -  A  one -dimensional 

quaal- steady  compressible  flow  model  is  used  to 
describe  the  gas  exchange  between  che  engine 


Fl  dM, 


dt 


9  dM 
dt 


(10) 


where , 


che  discharge  coefficient 

the  effective  flow  area 

the  upstream  stagnation  pressure 

che  downstream  stagnation  pressure  at 

che  effective  flow  area 

che  upstream  gas  temperature 

che  ratio  of  specific  heats. 


Due  to  che  large  pressure  difference  between 
Che  cylinder  and  che  exhaust  manifold  when  the 
exhaust  valve  just  begins  to  open,  the  flow 
through  the  exhaust  valve  can  become  choked.  In 
that  case,  when  sonic  conditions  prevail  at  che 
exhaust  valve,  the  critical  pressure  ratio  is 
evaluated  from: 


(12) 


COMBUSTIOH  PROCESS  •  The  combustion  model 
used  in  the  present  simulation  program  is  that 
developed  by  Uatson,  Piley,  and  Marzouk  [12],  It 
is  based  on  algebraic  expressions  describing  che 
fuel  bumig  rate  as  a  function  of  the  dominant 
controlling  parameters,  such  as  ignition  delay 
and  equivalence  ratio.  Chemical  reactions  are 
considered  to  commence  at  Che  dynamic  Injection 
point  and  to  cake  place  in  two  distinct  periods: 
an  ignition  delay  followed  by  a  period  of  heat 
release.  The  ignition  delay  is  considered  to  be 
che  time  interval  between  che  dynamic  injection 
point  (when  injector  needle  starts  to  lift)  and 
che  ignition  point  (start  of  heat  release  due  to 
combustion),.  The  total  length  of  che  ignition 
delay  is  related  to  che  cylinder  gas  temperature 
and  pressure,  the  engine  speed,  and  che  ignition 
quality  of  che  fuel  according  to  an  empirical 
expression  by  Hardenberg  and  Hase  [30]: 


IDa  -  (0.36+0,22Vp)  exp 


EA 


E  T_ 


— ) 

17190/ 


21.2 


P,  -  12.4 


0.63 


(13) 
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where, 

ID.  - 


the  ignition  da lay  in  crank  angle 

degrees 

cha  naan  piston  speed 

Cha  apparant  activation  energy 

618, 840/ (  Cn  ♦  25  ) 

fuel  catena  nuabar 

Cha  abaoluta  cylinder  pressure 

cha  absolute  cylinder  temperature . 


The  ignition  delay  in  seconds  can  be 
determined  by  simply  dividing  the  result  from 
Eq. (13)  by  the  quantity  6N.  where  N  is  the  engine 
speed  in  rev/nin. 

The  combustion  correlation  consists  of  two 
algebraic  expressions,  one  for  pre-aixed  fuel 
burning  and  the  ocher  for  diffusion  controlled 
fuel  burning.  •  The  total  instantaneous  fuel 
burning  rate  is  given  by  the  sum  to  the  two 
components; 


«t  *  Hp  +  Mi  (14) 

where  the  subscripts  t,  p,  and  d  denote  tocal, 
pre-aixed,  and  diffusion,  respectively.  In  order 
to  quantify  the  proportion  of  che  fuel  burnt  by 
either  mechanism,  a  'mode  of  burning  factor  ••  V 
is  used  to  express  che  cumulative  fuel  burnt  by 
pre-mlxed  burning  as  a  fraction  of  cha  cccal  fuel 
injected. 

a*  -  Mp  /  Mt  (15) 

The  mode  of  burning  factor  is  calculated  from: 

B,  -  1.0  -  0.746  •tr0-35  ID,*0  35  (16) 

where,  bcr  la  che  equivalence  ratio  of  che 
crapped  gas,  and  ID,  is  che  ign  icion  delay  in 
seconds.  Substituting  Eq.(lS)  into  Eq.(16),  che 
non-dimensional  fuel  burning  rate  can  be 
expressed  as; 


Mt  -  8,^(0  ♦  (1-B,,)  Md(t)  (17) 

For  che  pre-aixed  phase,  che  burning  race  is 
given  by: 


Mp(C) 


Cpl‘l 

Cpt  Cp2  T 


CP2*1 


(18) 


For  Che  diffusion  phase,  che  burning  race  is 
given  by: 


cd2'l  ':i2 

Kd(C)  -  Cdl  Cd2  T  exp(-Cdl  T  (19) 


where . 


£p2 

cdi 


-  2.0  +  0.002703  x  ID.2  -° 

-  5000 

-  10.0  *er 


•1.505 


-  0.48  CA{  0.«3 


The  variables  Cp2,  C_2,  Cdl.  and  Cd2  are 
experimentally  determined  'shape  factors'  which 
describe  Che  fuel  burning  rata  as  a  function  of 
engine  operating  parameters  [12], 


GAS -TO -WALL  HEAT  TRANSFER  IB  THE  CTLINDER- 
Kuaerous  correlations  for  the  gas-co-wall  heat 
transfer  coefficient  have  been  proposed  over  the 
years  (14],  In  the  present  model,  che  correlation 
by  Hohenberg  [39]  was  used.  Hohenberg  presented  a 
simplified  correlation  for  heat  transfer 
coefficient  which  made  allowance  for  the  rise  in 
gas  velocity  with  increasing  engine  speed  as  well 
as  che  variation  in  velocity  with  time  (crank 
angle)  expressed  as; 

hc  -  Cj.  Vd-0-06  p0.8  t-0.4  (Vp  +  c2)0-8  (20) 

where  che  constants  and  C2  are  empirical 
constants  which  can  be  adjusted  for  local 
variations  due  to  intake  swirl  or  combustion 
chamber  geometry.  Based  on  comprehensive 
measurements  on  various  types  of  DI  diesel 
engines,  Hohenberg  recommends  values  of  130  for 
Ci  and  1.4  for  C2.  Although  the  Hohenberg 
correlation  does  not  explicitly  Include  a  term 
for  radiative  heat  losses  in  che  cylinder,  the 
empirical  basis  of  Che  multiplying  constant 
Implicitly  allows  for  radiation. 

HEAT  TRANSFER  DDR  IK  INTAKE  AID  EXHAUST - 
The  heat  transfer  in  the  intake  and  exhaust  ports 
is  modeled  using  a  quasi-steady  heat  transfer 
correlation,  in  the  form  of  a  Nus suit -Reynolds 
number  relationship  given  by  Malchow  et  al  [15); 

Nud  -  K  ReD  a  (21) 

where  Raj)  denotes  the  Reynolds  number,  Uh 
subscript  D  denotes  the  diameter  of  the  Intake  or 
exhaust  pipe,  and  K  and  a  denote  empirical 
constants  taken  to  be  0.0774  and  0.769, 
respectively. 

CALCULATION  OF  METAL  BALL  TEMPERATURES  A 
simple  one -dimensional  heat  conduction  model  i  s 
used  to  evaluate  che  mean  inner  metal  surface 
temperature  of  che  combustion  chamber. 
Additional  heating  of  che  cylinder  wall  by 
friction,  and  cooling  of  che  wall  by  che  flow  of 
oil  over  che  surface,  is  neglected.  It  is  known 
from  experiments  chat  che  wall  cemperacure 
changes  in  the  actual  engine  are  very  small 
compared  Co  che  changes  in  gas  cemperacure.  The 
model,  therefore.  Cakes  this  fact  into  account  by 
updating  che  wall  Cemperacure  only  sc  the  end  of 
che  cycle  rather  Chan  at  each  step  during  che 
cycle. 

MECHANICAL  FRICTION  -  A  number  of  empirical 
correlations  have  been  proposed  to  model  total  or 
cycle  averaged  mechanical  friction  in  internal, 
combustion  engines  [16-23],  However,  there  are 
very  few  correlations  available  in  the  literature 
which  model  individual  frictional  components  on  a 
crank  angle  basis  rather  chan  on  a  cycle  basis. 
Rezeka  and  Heneln  [17],  developed  a  method  co 
determine  che  instantaneous  engine  frlccion 
during  firing.  The  method,  known  as  che  p-w 
method,  is  based  on  the  fact  Chat  che 
instantaneous  cylinder  gas  forces  and  the 
Instantaneous  frictional  forces,  inertia  and  1o.r<d 
forces  cause  che  instantaneous  variation  in  tUe 
flywheel  angular  velocity.  They  developed  a  sec 
of  correlations  to  calculate  che  lnstantem  ous 
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friction  for  each  of  cha  major  frictional 
components  which  were  ciaasiflad  under  two 
categories;  l)  Piaton  aaaaably  friction,  and  2) 
Crankcaae  aaaaably  friction. 

Platan  Assembly  Friction  -  Thraa  coaponanta 
of  the  piaton  aaaaably  friction  arc  considered: 
l)  Ring  Viscous  Lubrication  Friction  (Frvl)>  ^) 
Ring  Hixad  Lubrication  (Fr^) ,  and  3)  Piston 
Skirt  (Fpg): 

FRVl  “  »l[MVpwo(PgM+P#)]  -5(no+0. 41^)0  r  Tk  (22) 
fRML  *  e2irDnewc(PgM+Pe)(!-sin8)r  Tk  (23) 

Fps  “  *3(pVQLpr  T^/ho)  (24) 

Crankcase  Assembly  Friction  -  Crankcase 
friction  Include  three  aajor  components;  I)  Valve 
Train  Friction  (FyA^) ,  2)  Auxiliaries  and  Unload¬ 
ed  Bearing  Friction  (FAug),  and  3)  Loaded  Bearing 
Friction  (F^g): 


pVgL  “  GL,r  Tk  /Vw  (25) 


Faob  “  »5  #i  w 


(26) 


FLB  -  *6  -a-  D2  rc  P»«s  cose/Vw 
4 


(27) 


The  above  correlations  were  used  to 
calculate  the  instantaneous  friction  torque  (IFT) 
for  each  component.  The  total  instantaneous 
friction  at  any  crank  angle  is  given  by; 


IFT  -  2  IFTt  -  at  x  wt  [NM]  (28) 

where,  denotes  the  empirical  coefficients  for 
the  icb  friction  components,  and  w^  denotes  the 
functions  given  in  Eqs.(22)  through  (27).  The 
values  of  a^  and  w^  are  given  in  Table  1.  Figure 
3  shows  the  torques  for  Che  different 
instantaneous  frictional  components  calculated 
using  the  correlations  in  Table  1.  For  this 
figure,  the  engine  is  running  at  900  rpm,  no 
load.  The  summation  of  the  instantaneous  values 
of  all  frictional  com-  ponents  gives  the 
instantaneous  total  friction  shown  in  Figure  4. 

BLOtfBT  GAS  FLOW  •  The  blowby  model  used  in 
this  study  represents  the  flow  path  by  a  series 
of  chambers  (volumes)  connected  by  square-edged 
orifices  (ring  gaps)  as  shown  in  Figure  5.  The 
ring  gaps  are  assumed  to  be  the  only  openings 
for  the  gas  to  leak  out.  The  mass  flow  between 
these  volumes  can  be  expressed  in  terms  of  the 
pressure  ratio  and  the  stagnation  conditions.  The 
formula  is  for  one  -  dimens  iona  l .  unsteady, 
compressible  adiabatic  flow  through  an  orifice, 
wlch  a  constant  discharge  coefficient 

Mb  -  Cd  A,  Pi  TZ  gc  /  Rl  Tl  (29) 


where. 


ti 

2/y 

(i+D/i 

w 

■i )  \ 

*i  / 

—  Sing  viscous  lubrication 
-•*  Sing  mixed  lubrication 
-♦-Piston  skirt 

a 


Figure  3  -  Instantaneous  Friction  Components 
From  Hodel  Correlations  (Ref. 18) 
900  rpm.  No  Load 


Valve  train 

Auxiliaries  A 
unloadad  baarlngs 

— Loadad  bearings 


Figure  4  -  Total  Friction  Torqua  (Ref. 18) 
900  rpm,  No  Load 


C  d  la  the  discharge  coefficient  for  a 
square-edged  orifice  (approximately  equal  to 
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rotating  inertia  cor  qua ,  Troc,  la  calculatad  by; 

Troc  -  Xe  —  (32) 

dc 

where, 

I  -  che  moment  of  Inercla  calculated  from 
che  engine  geometry  and  masses  of  Che 
rotating  parts  (flywheel,  crankshaft, 
and  part  of  the  connecting  rod 
du/it  -  che  instantaneous  angular  acceleration 
1  d m2 


2  dd 
and  for  small  0, 


actual  sms  secs 

Figure  S  •  Modal  for  Blowby  Gas  Flow 
Past  che  Piston  Rings 

0.65).  The  mass  flow  (M^)  is  found  by  examining 
each  pair  of  adjacent  volumes  in  turn.  The  mass 
and  pressure  changes  within  each  volume  are 
related  by  che  perfect  gas  law: 


The  reciprocating  torque  (Trac) 
calculated  according  to: 


dTk 

Tk  w2 

— 

.  )  - 

dC 

d© 

dTk  . 

- - 

+ 

w2 -  ) 

(34) 

dt 

d© 

dMt,  - 


where , 

Vj  -  che  crevice  volume 

Twav  “  che  average  wall  temperature. 
Equations  (29)  and  (30)  can  be  numerically 
integrated  to  obtain  the  pressure  history  of  each 
volume,  and  more  importantly,  the  mass  flow  out 
of  the  main  cylinder  and  into  che  crankcase 
(blowby) . 

ENGINE  DYNAMICS  -  The  model  calculates  che 
instantaneous  angular  velocity  and  acceleration 
of  the  engine's  flywheel  during  cranking  and 
starting.  The  instantaneous  angular  velocity  of 
the  flywheel  depends  on  che  moment  of  inertia  of 
che  moving  parts,  and  on  che  ins cancaneous' 
difference  between  che  driving  torques  produced 
by  che  gas  pressure  and  the  scarcer  raocor  (during 
cranking) ,  and  the  resisting  torques  which  are 
Che  load  torque,  Che  inercla  torque,  and  the 
friction  torque.  Mathemac ic a  1  ly  ,  this 
relationship  can  be  expressed  as : 


”  Tscart  +  Tgas 


^fric  '  1 1 oad 


where , 

T^n,r  •  the  instantaneous  inertia  corque 
T,tart  “  ch*  instantaneous  scarcer  raocor 
corque 

T_M  -  che  instantaneous  gas  pressure  corque 
T®ric  -  che  instantaneous  friccion  corque 
Tioad  ”  the  Instantaneous  load  corque 

The  instantaneous  inertia  corque  (Tin#r) 
consists  of  two  part?;  1)  rotating  inertia 
torqua,  and  2)  reciprocating  inertia  corque.  The 


where , 

m  -  che  mass  of  the  reciprocating  parts 

(piston  assembly,  part  of  che  connecting 
rod,  etc.) 

r  -  che  crank  radius 
u  -  che  angular  velocity 
Tk  -  the  transformer  ratio 

sin©  r/1  sin©  cos© 

-  -  (35) 

[  1  •  (r/1  sin©)2  ]  5 


The  output  corque  produced  by  the  starter 
raocor  is  a  function  of  che  input  power  (battery 
voltage  and  current),  the  load  that  it  drives, 
and  che  mechanical  efficiency  of  the  starter 
These  determine  che  resulting  starter  speed.  The 
instantaneous  corque  output  of  the  scarcer  moior 
is  approximated  by  using  an  empirical  relation¬ 
ship  suggested  by  Poublon  et  al.  f 28 ]  as: 


Tstarc  “  Cs  **P(-Ks  N  > 


N  -  che  instantaneous  flywheel  speed  (in 
rev/min) 

CS,KS  -  are  empirical  constants  based  on  che 
scarcer  characteristics. 

The  value  cf  C,  and  K,  are  determined  by  plotting 
che  natural  log  of  the  s carter  motor  corque 
during  cranking  against  Che  starter  speed  for  a 
number  of  speeds.  The  resulting  points  are  curve 
fitted  using  a  least  squares  technique.  Then,  the 
ln(T,tart)  ls  ch,n  converted  back  to  its  linear 
form  and  che  constants  are  calculatad.  Figure  6 
shows  a  plot  of  the  scarcer  corque  calculated 
using  Eq. (36) . 
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Figure  6  -  Starter  Motor  Torque  aa  a 
Function  of  Engine  Speed 

The  instantaneous  gas  torque  (T-4J)  is 
computed  at  each  crank  angle  from  the  gas 
pressure  (P»ms) ,  piston  area  (Ap^s) ,  crank  radius 
(r),  and  length  of  the  connecting  rod  (1) 
according  to : 

Tgas  “  rTk(pgas  Apis  *  “rac  8  sine)  (37) 

vhere, 

®r,c  "  the  of  the  reciprocating  piston 

assembly 

Tk  -  the  transformer  ratio. 

The  Instantaneous  load  torque  (T^oad)  can  be 
determined  from  the  characteristics  of  the  load. 
Figure  7  shows  a  plot  of  the  instantaneous  torque 
components  during  steady  state  motoring, 
calculated  using  Eqs.(31)  through  (37). 
Substituting  Eqs.(32)  and  (34)  into  Eq.  (31)  and 
re-arranglng  the  resulting  expression,  the 
angular  acceleration  of  the  engine  can  be 
calculated  by: 


CRANK  ANCLE 


^start+'pgae'^fric'^load'l*r^"^k  d u 
dt  mr2  Tk2  +  Ie 


By  substituting  Eq.(38)  into  Eq.(33),  and 
re-arranging  the  resulting  expression,  the 
instantaneous  angular  velocity  of  the  flywheel  at 
the  next  crank  angle  (time)  step,  can  be 
approximated  by; 


2  2 
"1+1  •  "1 

dt 

2  AO 

or. 

r  **  s 

"1+1  “  2  AS  -  +  "I'M 

(39) 

L  dc  J 

Figure  8  shows  the  instantaneous  angular  velocity 
and  acceleration  calculated  from  Eqs.(39)  and 
(38)  for  a  single  cylinder  diesel  engine  running 
at  1100  rpm,  no  load. 


Figure  8  -  Instantaneous  Angular  Velocity  and 
Acceleration,  1100  rpm.  No  Load 


Figure  7  -  Instantaneous  Torque  Components 
1100  rpm,  No  Load 
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METHOD  OF  SOLOTICM 

For  Crons  lone  operation,  ehs  differential 
•qusclons  for  etas  ttarss  control  volumes  oust  bs 
solvsd  ss  sn  initial  vslua  problem  sines  ehs 
conditions  se  ehs  snd  of  ehs  cycle  will  not 
( necessarily)  equal  ehose  se  the  beginning.  The 
decsils  of  the  solution  are  given  in  reference 
[11].  The  starter  motor  torque  is  applied  eo  the 
engine  as  driving  function  which  produces  torque 
labslances  in  the  systea  and  changes  in  the 
systea's  operating  conditions.  These  changes  are 
computed  from  the  elae  chat  the  starter  key  is 
turned  on.  until  the  engine  fires  and  accelerates 
to  lea  idling  speed,  or  until  a  specified  cine 
(nuaber  of  cycles)  has  expired.  The  engine  is 
assumed  to  have  started  when  Che  engine  speed 
reaches  or  exceeds  a  specified  rpa.  The 
calculation  is  then  stopped,  and  the  results  are 
printed. 

The  predictor-corrector  aethod  used  in  the 
present  study  consists  of  Che  following  formulas 
[40]: 

Predicting  equation: 

Yi+l  -  ♦  h  f'(X1,Yl)  (40) 

Correcting  equaelon: 

Y'i+1  *  Y1  ♦  h/2[f‘(Xi.Y1)+f'(X1+1.Y1+1)]  (41) 
where, 

Yj  -  the  value  at  the  current  step 
Y^+i  -  the  predicted  value  at  the  next  step 

Y’  -  the  corrected  value  at  the  next  step 

h  -  the  step  size. 

EXPERIMENTAL  APPARATUS  AMD  PROCEDURES 

The  engine  used  to  obtain  experimental  data 
for  comparison  with  the  model  predictions  was  a 
F1L  210  DEUT2  MAG  single  cylinder,  four-stroke, 
air-cooled,  direcc  injection  diesel  engine.  The 
engine's  specifications  are  given  the  appendix. 
The  engine  was  equipped  with  a  starting  motor  and 
battery  for  starting  the  engine  without  a 
dynomoaeter.  and  for  motoring  the  engine 
(i.e.  cranking  without  combustion). 

A  series  of  starting  tests  were  conducted 
with  the  engine  warmed-up  to  idling  temperatures. 
The  inlet  air  temperature  and  pressure  were  88  °F 
and  29.7  in  Hg. ,  respectively  The  oil 
temperature  was  130  °F.  Fuelling  was  0142  mm  per 
cycle  (F/A  -  0.02).  Engine  behevlor  was  monitored 
from  the  instant  the  scarcer  key  was  turned  on, 
referred  to  as  the  starter-on  point,  to  the  time 
when  the  engine  started  and  accelerated  to  reach 
its  final  idle  speed  of  approximately  l 100  rpm. 

The  cylinder  gas  pressure  was  measured  using 
an  AVL  12QP  300  CVK  pressure  transducer,  flush 
mounted  in  the  cylinder  head.  The  transducer 
signal  was  fed  into  a  Klstler  modal  304  E  charge 
amplifier,  and  its  output  displayed  on  an 
oscilloscope.  The  crank  angle  marks  used  for 
referencing  Che  pressure  were  measured  and 
recorded  simultaneously  using  a  magnetic  pick-up 
transducer  mounted  opposite  to  the  93 -tooth  ring 


gear  on  £he  angina's  flywheel.  A  separate 
magnetic  pick-up  was  used  to  indicate  TDC. 

The  angular  velocity  of  the  flywheel  was 
computed  by  measuring  the  time  taken  to  move 
through  a  fixed  angle  of  rotation.  An 
electro-magnetic  pick-up  was  installed  in  the 
flywheel  housing  opposite  the  teeth  in  such  a  way 
to  produce  a  pulse  signal  to  indicate  the 
presence  of  a  tooth  on  the  flywheel  ring 
gear.  The  instantaneous  angular  velocity  at  the 
flywheel  end  was  calculated  from: 

ii>£  «  AS  /  tk  (42) 

where . 

4©  -  the  angle  between  two  successive  gear 
teeth 

-  2*  /  Mt 

Nc  -  che  number  of  teeth 

tk  -  the  time  taken  to  move  4© 

-  Yk  /  fck 

Yfc  -  the  number  of  clock  counts  during  40 
fCk  “  the  clock  frequency. 

Substituting  che  expressions  for  40  and  tk  into 
Eq . (42)  gives  che  following  expression  for  the 
angular  velocity: 


«f  -  2*  fclt  /  Yk  Nt  (43) 

The  power  input  to  the  starter  motor  vtu 
determined  by  measuring  the  battery  voltage  and 
Che  starter  current  while  cranking  the 
engine.  The  input  power  was  than  calculated  by; 

Ps  -  V,  x  I,  (44) 

where, 

Ps  -  che  input  power  to  drive  the  starter 
in  watts 

Vs  -  the  bacery  voltage 

Is  -  che  scarcer  current. 

COMPARISON  BETWEEN  EXPERIMENTAL  AND  MODEL  RESULTS 

Figure  9  shows  a  comparison  of  che 
analytical  and  experimental  instantaneous  angular 
velocity  and  cylinder  pressure  as  a  function  of 
crank  angle  for  che  first  12  cycles  during  the 
starting  process.  Initially,  the  model  predicts  a 
slightly  higher  angular  velocity  chan  the 
experiment.  The  lower  experimental  angular 
velocity  is  probably  due  to  mechanical  losses  in 
the  starter  which  were  not  accounted  for  in  the 
model.  The  starter  off  speed  for  this  test  was 
about  90  rad/sec  (8S9  rpm).  Once  the  starter 
torque  vas  turned  off  in  che  model,  che 
calculated  angular  velocity  approached  che 
measured  angular  velocity.  By  Che  CIsm  the  engine 
had  accelerated  to  its  idle  speed  of  about  120 
rad/sec  (1150  rpm),  Che  angular  velocities  from 
che  cast  and  from  che  model  agreed  very  closely. 

In  che  first  three  firing  cycles,  the 
calculated  peak  pressurss  were  significantly 
higher  than  Che  measured  pressures.  This  may  be 
partly  due  to  che  assumptions  used  in  the 
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Figure  9  •  Instantaneous  Angular  Velocity  and 
Cylinder  Pressure  During  Cranking  and 
Starting  Iron  the  Model  and  Experiment 


combustion  nodal.  Once  the  nodal  determines  that 
conditions  are  favorable  for  Ignition  to  occur; 
chat  Is,  the  ignition  delay  is  shorter  chan  Che 
available  Cine,  It  asaunes  complete  coabustlon  of 
all  Che  fuel  Injected.  In  Che  actual  engine, 
however,  this  is  not  the  case.  Partial 
coabustlon,  and  sonatinas  no  coabuscion  (aisflre) 
can  occur,  particularly  under  marginal  starting 
conditions.  Hence,  the  higher  pressures  predicted 
by  the  nodal  nay  have  resulted  from  the  higher 
.energy  release  froa  Che  additional  fuel.  The 
higher  cylinder  pressures  would  also  account  for 
the  higher  angular  velocities. 

There  appears  to  be  a  phase  shift  in  time 
(crank  angle)  between  the  two  curves  This  may  be 
due  to  a  possible  error  In  measurement  of  TDC  in 
Che  engine,  or  a  slight  difference  in  the 
starting  crank  angle  used  In  the  comparison. 

In  susuiary,  the  comparisons  between  the 
analytical  and  experimental  results  appears 
satisfactory  for  the  starting  conditions 
Investigated.  Uhlle  there  is  some  difference  in 
the  absolute  values  predicted  by  the  model,  the 
relative  trends  seem  correct. 

MODEL  LUDTAXIOKS 

The  model  Is  necessarily  cailored  quite 
extensively  to  the  cast  engine,  both  in  the 
physical  constants  describing  the  engine  and  to 
the  dynamic  and  thermodynamic  models  used.  The 


empirical  basis  for  the  coabustlon  model,  starter 
model,  and  valve  flow  coefficients  limits 
usefulness  of  the  results  obtained  froa  the 
simulation  only  to  the  test  engine.  However,  the 
methodology  used  to  develop  the  model  can  be 
applied  to  other  engines  with  minimum  difficulty. 

Several  areas  in  Che  model  need  further 
development.  First  of  all,  Che  effects  of 
friction  heating,  and  oil  cooling  on  the  cylinder 
walls  during  cranking  and  starting  may  be 
significant,  and  should  be  accounted  for  in ‘the 
model.  Secondly,  improvement  Is  needed  In  starter 
model  to  include  load  effects  on  the  starter 
output  during  Initial  cranking.  Also,  the 
constants  used  in  the  scarcer  model  should  be 
adjusted  to  account  for  friction  In  the  starter. 
Thirdly,  extension  of  che  model  for  multi- 
cylinder  engines  would  be  useful  In  studying  che 
effect  of  several  cylinders  on  engine 
startabillcy. 

CONCLUSIONS 

The  comparisons  between  che  analytical  and 
experimental  results  generally  show  that  the 
model  satisfactorily  describes  che  starclng 
behavior  of  the  test  engine  at  room  temperatures. 
However,  cold  room  experiments  definitely  needed 
to  validate  model  predictions  for  startabillcy  at 
low  ambient  temperatures. 

The  basic  assumptions  employed  In  the 
'filling  &  emptying*  approach  to  cycle  analysis 
can  be  extended  to  the  low  rpm  range  to  study 
engine  behavior  during  initial  cranking  and 
starting. 

The  inclusion  of  detailed  calculations  of 
instantaneous  friction,  blowby,  and  heat  transfer 
coupled  with  engine  dynamic  calculations  in  a 
cycle  simulation  program  provides  an  Interesting 
framework  in  which  to  study  the  combined  effects 
of  changes  In  CR,  ambient  temperature,  and 
fuel/air  ratio,  etc  on  engine  startabillcy. 


- Description _ 

effective  flow  area 
area  of  piston 

total  surface  area  of  combustion 
chamber 

optimum  coefficients  for  friction 
model 

mode  of  burning  factor 
discharge  coefficient 
diffusion  shape  factors  for 
combustion  model 
cetane  number 
specific  heat 

pre-mlxed  shape  factors  for 
combustion  modal 

compression  rstio  / 

Starter  coefficients  \ 

diameter  of  cylinder  ^S, 

activation  energy 

battery  voltage 

friction  due  to  auxiliaries  & 

basrlngs 


NOMENCLATURE 


Svmbol 

Aeff 


al-a6 


3. 

=d 

cdl-cd2 


CP 

CP1’CP2 

CR 

Cs,Ks 

D 

EA 

Es 

faub 
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Fck 

FUS 

FMEP 

FPS 

fRML 

frvl 

Fs 

fval 

G 

5c 

h 

hc 

hf 


10. 

10, 

I. 

I  FT 

I. 

L 

5l 

«. 

«b 

«i 
* 


rec 


Mu 

N 

nC 

no 


l8*« 


Q 

R 

r 

rc 

Ih 

Tc 

TOC 

Tfrlc 

Tiner 

ck 


Tk 

Tload 

T. 

Trec 

*rot 

Tstart 


*w»v 

u 

vl 

V„ 


clock  frequency 

friction  duo  co  loaded  bearing, 
friction  naan  affacciva  prasaura 
friction  due  to  piston  skirt 
friction  due  co  ring  nixed 
lubrication 

friction  due  co  ring  viscous 
lubrication 

stoichioaatric  fuel/air  ratio  . 
friction  due  co  valve  train 
nuaber  of  valves  par  cylinder 
gravitational  constant 
specific  enthalpy 

cylinder  gas- to- veil  heat  transfer 

coefficient 

enthalpy  of  formation 

oil  file  thickness 

crank  angle  step  sire 

Ignition  delay  in  degree 

ignition  delay  In  seconds 

aoment  of  inertia  of  angina 

Instantaneous  friction  torque 

scarcer  current 

length  of  connecting  rod 

length  of  piston  skirt 

valve  spring  load 

■ass  of  air 

■ass  of  blowby  gas 

■ass  of  diffusion  burned  fuel 

■ass  of  fuel 

■as a  in  control  voluae  1 

■ass  of  pre-alxed  burned  fuel 

■ass  of  reciprocating  parts 

total  aass  of  engine  parts 

■ass  of  fuel  unburned 

engine  speed,  rpa 

nuober  of  coaprasslon  rings 

nuaber  of  oil  rings 

nuober  of  flywheel  gear  teeth 

Nusselt  nuaber 

elastic  ring  pressure  forces 

pressure  due  co  cylinder  gases 

power  input  to  starter 

heat  transfer  race 

universal  gas  constant 

crank  throw  radius  (1/2  stroke) 

radius  of  journal  bearing 

Reynolds  nuaber 

gas  constant 

revolutions  per  alnute 

cylinder  gas  teaperature 

cop  dead  canter 

total  friction  torque 

torqua  due  to  inertia 

Ciae  between  successive  flywheel 

gear  teeth 

transformer  ratio 

torque  due  co  load 

mean  teaperature 

torque  due  to  reciprocating  parts 
torque  due  to  rotating  parts 
torque  due  co  starter  aotor 
cylinder  wall  teaperature 
average  cylinder  wall  teaperature 
specific  Internal  energy 
voluae  of  control  voluae  1 
clearance  voluae 


*1*1 

yk 

x 

u 

“f 

* 

M 

u 

8 


cylinder  voluae  swept  by  piston 

velocity  of  the  gas 

piston  velocity 

Voltage  Input  co  Scarcer 

width  of  compression  ring 

width  of  oil  ring 

corrected  value  at  next  step,  1*1 

value  at  current  step  1 

predicted  value  at  next  step 

nuaber  of  clock  counts  per 

flywheel  tooth 

ratio  of  r/L 

angular  velocity 

angular  velocity  of  flywheel 

density  of  Inlet  air 

dynaaic  viscosity 

klneaatlc  viscosity 

crank  angle  In  degrees 

ratio  of  specific  heats  ( cp/c v) 

equivalence  ratio 

equivalence  ratio  of  trapped  gas 
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Engine  Specification*  for 
Dautx  FIL  210  MAC 


TabU  l 


Combustion  Systan .  D.I. 

No.  of  Cylinders  .  1 

Boro  x  Stroko  (an)  .  95  x  95 

Dlsplacanent  (ce)  .  673 

Compression  Ratio  (non.)  .  17:1 

Clearance  Volume  (ce)  .  0.4208 

Piston  Skirt  Langth  (an)  . 66.5 

Piston  Sklrt/Linor  Clearance  (on)  .  0.031 

Connaetlng  Rod  Langth  (an)  .  160 

Journal  Baarlng  Dlaaatar  (an)  .  95.19 

No.  of  Conpraaslon  Rings  .  2 

Width  of  1st  Ring  (am)  .  2.54 

Width  of  2nd  Ring  (on)  .  2.54 

No.  of  Oil  Rings  .  1 

Width  of  Oil  Ring  (an)  .  3.60 

Voluna  Batvaan  Rings  162  (cc) .  2.621 

Araa  Batvaan  Rings  1  &  2  (caj) .  0.0222 

Voluns  Batvaan  Rings  2  &  0  (cc)  .  1.818 

Araa  Batvaan  Ring  2  6  0  (  )  .  0.0222 

Mass  of  Piston  &  Rings  (kg)  .  1.18 

Mass  of  Connacting  Rod  (kg)  .  0.6S 

Rotating  Monant  of  Inertia  (kg-a-s)  ....  0.384 

Are*  of  Cylindar  Head  (cm2)  .  75.0 

Are*  of  Piston  (cn2)  .  99.15 

Araa  of  Linar  at  TDC  (cn2)  .  0.358 
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Injection  lining  (Dag.  BTDC)  .  23 

Intaka  Valve  Opens  (Dag.  BTDC)  .  21 

Intake  Valve  Close  (Dag.  ABDC)  .  62 

Exhaust  Valve  Opens  (Dag.  BBDC)  .  62 

Exhaust  Valve  Close  (Dag.  ATDC)  .  21 


SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 
6 . 1  Conclusions 

1.  The  basic  function  of  the  cetane  scale  is  to  rate  the 
autoignition  behavior  of  a  fuel  by  matching  its  ignition 
delay  to  that  of  a  blend  of  two  primary  reference  fuels 
under  specified  conditions  of  inlet  temperature,  engine 
speed,  coolant  temperature,  fuel-air-ratic,  injection  timing 
and  injector  opening  pressure.  The  only  engine  variable 
is  the  compression  ratio.  The  cetane  scale  was  not 
intended  to  rate  any  other  combustion,  performance  or 
emissions  parameter. 

2.  Three  intrinsic  problems  exist  in  the  cetane  rating  of 
fuels . 

a-  Inability  to  rate  low  ignition  quality  fuels.  The 
solution  is  to  raise  the  inlet  air  temperature 
above  the  150  F  specified  in  the  ASTM  D-613  proce¬ 
dure  . 

b-  Instability  of  the  ignition  delay  meter  in  rating  low 
ignition  quality,  fuels .  This  has  been  found  to  be 
due  to  the  cycle-to-cycle  variation  in  the  I.D.  of 
these  fuels  with  the  meter  trying  to  reflect  what 
is  actually  occuring  in  the  engine.  No  solution 
is  known  at  this  time  to  the  problem  of  cycle-to-cycle 
variation  in  the  I.D. 

c-  Inconformity  of  the  CFR  combustion  system  with  actual 
diesel  combustion  systems.  A  possible  solution  is 
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to  redesign  the  combustion  system.  Part  of  the  pro¬ 
blem  may  be  solved  for  the  existing  system  by  opti¬ 
mizing  the  injection  pressure  such  that  the  rating 
of  the  fuel  is  independent  fo  the  needle  opening 
pressure. 

Three  fuel/ engine  performance  problems  are  examined: 
a^-  Invalidity  of  C.N.  in  correlating  the  I.D.  for  dif¬ 
ferent  fuels  in  different  engines  under  all  operating 
conditions.  The  I.D.  should  be  correlated  in  terms 
of  three  basic  parameters:  mean  temperature,  mean 
pressure,  and  a  global  activation  energy.  The  glohal 
activation  energy  can  be  correlated  with  the  C.N. 
b-  Inaccuracy  of  C.N.  in  rating  the  overall  performance 
and  emissions  behavior  of  fuels  in  diesel  engines. 
This  correlation  should  be  done  in  terms  of  I.D. 
and  the  combustion  parameters  which  characterize 
the  premixed  and  diffusion  combustion  regimes, 
c-  Inability  of  C.N.  in  rating  the  cold  startability 
of  diesel  engines  on  different  fuels.  Startability 
depends  on  the  success  of  both  the  autoignition  and 
the  subsquent  combustion  reactions.  Many  parameters 
other  than  C.N.  need  to  be  included. 

The  C.N.  alone  cannot  cate  the  cold  startability  of 
diesel  engines  on  alternate  fuels.  Other  fuel  properties 
play  a  major  role  in  the  starting  of  engine.  The  most 
Important  property  ts  fuel  volatility.  This  conclusion 


is  supported  by  steady  state  tests  on  the  misfiring 
temperatures  and  ignition  delay  periods.  Also  it  is 
supported  by  actual  starting  tests  at  low  ambient 
temperatures  in  a  cold  room. 

6 . 2  Recommendations 

1-  Study  the  cetane  rating  at  inlet  air  temperatures  higher 
than  the  150  P  specified  by  the  ASTM  D-613  procedure. 
Determine  the  most  suitable  temperature  whi±h  allows 

a  suitable  size  of  electric  heater  to  be  used. 

2-  Investigate  the  phenomena  of  cycle-to-cycle  variation 
in  the  autoignition  process  of  fuels. 

3-  Investigate  the  effect  of  optimizing  the  needle  opening 
pressure  for  each  fuel  on  the  cetane  rating. 

4- r  Develop  correlations  between  the  C.N.  and  the  global 

activation  energy  for  the  preignition  reactions  for 
different  fuels,  to  be  used  for  correlating  I.D.  with 
C.N.,  temperature  and  pressure. 

5-  Develop  correlations  between  I.D.  and  the  parameters 
for  the  premixed  and  diffusion  combustion  and  the  fuel/ 
engine  performance  parameters.  Relate  these  parameters 
to  C.N.  by  using  C.N. -I.D.  correlations. 

6-  Make  a  combustion  failure  analysis  during  engine  starting 
on  different  fuels  and  ambient  temperatures.  This  will 
Identify  the  engine  and  fuel  parameters,  other  than 

C.N.,  which  affect  engine  startability . 

\ 

7-  Investigate  this  process  of  spray  formation,  autoignition 
and  combustion  at  low  ambient  temperatures. 
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